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ABSTRACT "E 


The innovative idea of distributing the functionality of current larger satellites 
among smaller, cooperative satellites has been sincerely considered for assorted space 
missions to accomplish goals that are not possible or very difficult to do with a single 
satellite. Additionally, the utilization of smaller satellites 1s maximized within 
formations and clusters to conduct missions such as interferometry and earth-sensing. 
This paper presents a methodology to describe, populate and analyze numerous 
formation designs employing the use of Hill's equations of motion to describe a 
formation’s dynamics. These equations of motion are then programmed into a 
MATLAB code to produce Cartesian elements for input into a Satellite Tool Kit" 
(STK) simulation that demonstrates numerous possible cluster formation designs. 
These simulations are then used to determine AV requirements for overcoming LEO- 
type perturbations that were modeled within STK's High Precision Orbit Propagator 
(HPOP). 


Finally, components from two subsystems [Attitude Determination and Control 
(ADCS) and Propulsion], using the AV calculations from the simulation analysis and 
current advances in MicroElectroMechanical systems (MEMs) and nanosatellite 


technology, are presented based on a mass constraint of 10kg for the entire satellite. 
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I. INTRODUCTION 


Nanotechnology provides the capability to manipulate matter at the atomic 
level. In the future, we will measure the way we design and build our 
systems by the atom, not by the pound. Today, we are developing material 
systems, at the molecular level, that are 100 times stronger than steel at 
1/6 the weight. We will also develop sensors and detectors capable of 
responding to a single photon of light or the stimulus from a single 
electron. Using nanotechnology, we will build systems on a scale 1000 
times smaller than today -- at true molecular level. They will be based on 
concepts emerging from biology, quantum mechanics and chemistry, all of 
which have no current parallel. [Ref 1] 


A. NANOTECHNOLOGY PHENOMENON 


Commercially and militanly, space systems based on MicroElectroMechanical 
systems (MEMS), nanoscale design and matenals, low power quantum electronics, and 
high bandwidth photonics are of special interest, as are the demonstrations of space 
subsystems based on these technologies. Significant reductions in individual spacecraft 
mass and cost, large gains in capabilities and robustness, and novel architectures 
involving large constellations and closely coupled spacecraft are expected with the 
introduction of these technologies. Novel spacecraft architecture concepts for Earth 
orbiting missions including communication, navigation, remote earth-sensing, and 
monitoring of the local space environment that can be created with the introduction of 


micro/nano-technologies will be outlined briefly below. 


B. DEFINING THE STARTING LINE 


First of all, it is worth defining what we mean by a small satellite. The spirit of the 
current small satellite world is encompassed by the slogan ‘Faster, Better, Smaller, 
Cheaper’. Small satellite projects are characterized by rapid development when 
compared with the conventional space industry, often ranging from six to thirty-six 


months. Leading-edge technology is routinely included in order to provide innovative 


solutions, permitting lighter satellite systems to be designed inside smaller volumes. 
Frequently, traditional procedures, with roots in the military and manned space programs, 
can no longer be justified, and low cost solutions are favored to match the reducing space 
budgets. So in many ways it is the philosophy, and not the size or mass of the satellite 


that matters. 


Many terms are used to describe this rediscovered class of satellites, including 
SmallSat, Cheapsat, MicroSat, MiniSat, NanoSat and even PicoSat. The US Defense 
Advanced Research Projects Agency (DARPA) refers to these as LightSats, the U.S. 
Naval Space Command as SPINSat's (Single Purpose Inexpensive Satellite Systems), and 
the U.S. Air Force as TACSat's (Tactical Satellites). Nevertheless, in recent years a 
general method of classifying satellites in terms of deployed mass has been generally 


adopted (see Table 1.1). 


The boundanes of these classes are an indication of where launcher or cost 
tradeoffs are typically made, which is also why the mass is defined including fuel (i.e. 
"wet mass”). Although the satellites in the 500-1000kg are typically designated as a ‘small 


satellite’, this causes confusion and until a better term appears it will be defined here as a 


medium sized satellite. 
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Table 1.1 Satellite Classifications 










Small Satellites 





The mass distribution for small satellites plotted in Fig. 1.1 illustrates that there 
re no clear mass boundaries, although there is a general lack of spacecraft in the 100- 
'DOkg class. The positive sloping line in the figure depicts a gradual increase in the mass 
‚f launched satellites. Although this trend will continue as heavier satellites are launched 
vithin the next few years, the number of smaller satellites launched should start to 


emove this trend. 
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Figure 1.1 Small Satellite Launch Mass 


E. CURRENT AND FUTURE EMPHASIS 


Few modern Nano and Pico-Satellites weighing less than 10kg have been 
aunched, although there is considerable interest in this area as advanced microsat 
echnology is applied to miniaturize satellite systems even further. Some launched 


nicrosats almost fall into this category with weights of 11-14kg, most notably the 


AMSAT microsat-series [Ref 2]. These satellites are cubical in shape and measure less 


than 150mm on each side. 


Nanosatellites are attractive to many educational institutions to get involved in 
space, as commonly available technology now makes this type of satellite feasible and 
most importantly affordable. Although operational picosatellites weighing less than Ikg 
are still some time off, Aerospace Corporation just recently launched a set of 
picosatellites (see Chapter 2A: JAWSAT below). Increasingly, micro and 
nanotechnology makes it possible to fabricate entire satellite sub-systems, and possibly 
entire satellites on a chip. Considerable effort is being spent on these Femtosatellites 
weighing less than 0.1kg, with applications in remote inspection, distributed 


measurement and disposable sensors. 


For Nanosatellites, autonomous operation using a single on-board computer is 
feasible, making use of technology developed for laptop and. palmtop computers. To 
minimize mass, active attitude and orbit control are often ignored, and omni-directional 
antennas are employed. The main limits are set by the downlink and power generation 
systems. The downlink data rate is limited by the orbit average power generation, and has 


to be operated at low data rates, or in burst mode. 


The key to the successful development of nanosatellites and constellations is 
advanced technology. Each of the many enabling technologies toward nanotechnology 
represents a breakthrough in performance, capability, or application in a unique way. 
This technological challenge is formidable, since currently, the smallest "full-service" 
microspacecraft weighs 100 kg (220 Ib) or more - size must be reduced by almost a factor 
of five. This will require revolutionary advances in microelectronics and spacecraft 


component technologies. 


In addition, at present day, almost all space missions are flown as single 
spacecraft. This 1s because controlling spacecraft in flight 1s a very complex process. The 
problem of flying several spacecraft as one system 1s further compounded by the complex 
communication path from a constellation of spacecraft in flight high above Earth with 
communication stations on the ground. In order to work properly, the spacecraft will have 


to behave “intelligently” -- autonomously staying in constant contact with each other, 
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sharing information, and re-configuring onboard instruments and systems to behave as a 


single unit. 


D. SIZE OBSTACLES 


A fundamental problem in spaceflight is spacecraft size. It is very expensive to 
place satellites into space. For example, using expendable cost with eight launches a year 
for the Space Shuttle, it costs about $13,200 per kilogram ($6,000 per pound) to deploy a 
payload into low-Earth orbit. As a result, engineers try to design spacecraft to be as small 
as possible. It has been estimated that satellites using nanotechnology could measure 15 
inches wide by 2 inches thick and weigh about 2 pounds. These satellites would then 


require smaller and therefore less expensive launchers. [Ref. 3] 


Most Earth-observing spacecraft with science payloads weigh in excess of 1,000 
kg (2,200 lbs). Microsatellites are much smaller, typically in the 100 kg (220 Ib) or larger 
range. Nanosatellite are even smaller -- in the sub-20 kg (44 lbs) range. Small spacecraft 
are nothing new. In fact, the first satellite launched into space, Sputnik, weighed only 
about 90 kg (180 lbs). However, small spacecraft even today are very limited in their 
capabilities. Typically they lack any propulsion, possess only limited ability for attitude 
control, and carry one single-function payload. 


What is needed is a new era of "smart", miniature spacecraft that will be "full 
service", meaning they will carry a wide range of spacecraft services including guidance, 
navigation and control, attitude control, propulsion, high bandwidth and complex 


communication functions. 


Nanosatellites not only reduce launch costs, but they also reduce the risks 
associated with flying missions. Currently, several instruments and payloads are flown on 
a single, large spacecraft. A single instrument or system failure may severely degrade or 
disable the entire mission. Constellations of numerous spacecraft, each carrying 
complementary instruments reduce the risk of an entire mission failing if one system or 


instrument fails. 


Е. CONCEPTS OF THE AEROSPACE CORPORATION 


Most of the radical new concepts for building and using spacecraft represented by 
the nanosatellite concept were developed at The Aerospace Corporation and formally 
introduced in a paper, "The Concept of 'Nanosatellite' for Revolutionary Low Cost Space 
Systems," presented at the 44th International Astronautical Federation Congress in Graz, 


Austria, in 1993. 


Authors included Robinson, Siegfried Janson, Ph.D., -- who coined the term 
nanosatellite and originated the concept [Ref 3] -- and Henry Helvajian, Ph.D., a senior 
scientist in the microtechnology center and editor of the just-published book, 


"Microengineering Aerospace Systems.” 


A series of reports written by Helvajian, Janson, and Robinson and issued after 
the 1993 conference presented the details on how to design, build, power and maneuver 
nanosatellites. These nanosatellite technologies are now being explored by a number of 
national and international research organizations in addition to The Aerospace 


Corporation 


The First International Conference on Integrated Micro/Nanotechnology for 
Space Applications was hosted jointly by the Aerospace Corporation and Johnson Space 
Center in Houston, from October 30" through November 3'* 1995. The purpose of the 
conference was to bring together scientists and engineers from the fields of 
microtechnology, nanoelectronics and space technology to explore the possibilities for 
applying newly emerging capabilities in microtechnology to space operations. The 
evolution of microelectronic technology coupled with the growth of MEMS in the past 4- 
5 years has had a significant impact in the commercial terrestrial sector. This influence 
can be evidenced particularly in sensor, optical switching and mass data storage 
applications that have been inserted into major industries such as transportation, 
medicine, telecommunications and computers. The focus of this conference was to 
anticipate and extend the incorporation of nanoelectronics and MEMS into Application 
Specific Integrated Microinstruments (ASIMs) in order to revolutionize the development 


of space systems. 


Е. THE SCOPE 


This paper sets out to introduce, explore and design nanosatellite formation 
designs with their required orbital dynamics. Chapter two will introduce past and current 
systems that are based on novel nanosatellite concepts. Numerous systems will be 
covered in areas ranging from experimental military designs to innovative student-driven 
ideas and future-looking commercial enterprises, all expecting to capture the quickly 


developing realm of ‘smaller, faster, cheaper’ with regards to nanosatellite technology. 


This paper presents a methodology to describe, populate and analyze numerous 
formation designs employing the use of Hill’s equations of motion to describe a 
formation’s dynamics. These equations of motion are then programmed into a MATLAB 
code to produce Cartesian elements for input into a Satellite Tool Kit" (STK) simulation 
that demonstrates numerous possible cluster formation designs. After utilizing 
MATLAB and STK to create formation simulations modeled within STK’s High 
Precision Orbit Propagator (HPOP), low-earth orbit (LEO) perturbations will be analyzed 
to understand the required delta-V (AV) to maintain the formation within given 


dimensional criteria. 


Finally, components from two subsystems [Attitude Determination and Control 
(ADCS) and Propulsion], using the AV calculations from the simulation analysis and 
current advances in MicroElectroMechanical systems (MEMs) and nanosatellite 


technology, are presented based on a mass constraint of 10kg for the entire satellite. 
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II. LITERATURE REVIEW: PAST AND CURRENT 
SYSTEMS 


A. PAST SYSTEMS 


1. 1958 


The first U.S. earth satellite, Explorer I, was almost the first nanosatellite by 
today’s standards. Launched by a modified Army Ballistic Missile Agency Jupiter-C on 
February 01, 1958 it, containing the upper stage, measured 2.03m long and 150mm in 
diameter, and with a mass of 13.6kg and an orbit of 347x1,859 km at 33.2 deg 
inclination. Explorer I, developed by the Jet Propulsion Laboratory, carned the U.S.-IGY 
(Intemational Geophysical Year) experiment of James A. Van Allen and resulted in the 
discovery of the radiation belt around the earth [Ref. 4]. Although most US satellites to 
follow would weigh more and become bigger in size, there was one that could be 


considered the ‘forefather’ of today’s nanosatellites: Vanguard-1. 


Although the first official ‘nanosatellite’ would not appear for another thirty 
years, the first satellite launched in the 1-10kg mass range is ironically the oldest satellite 
still orbiting Earth. Originally, a simple nose cone was to be carried on Vanguard-1, but 
in July of 1957 it was decided that a small 1.47-kg (3.25 pound) test satellite would be 


used instead to exercise the tracking stations (see Fig. 2.1). 
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Figure 2.1 Vanguard-1 


Launched March 17, 1958 into an elliptical (656x3866 km), 34.25 deg inclination 
orbit, this test satellite would consist of a simple 16 centimeter (6.4 inch) polished 
aluminum alloy sphere equipped with two transmitters operating at frequencies around 
108 MHz. Although the Vanguard 1 satellite must be looked at as a test satellite its 
technical and scientific results 1s impressive. Technically its purpose was to test апа 
evaluate solar cells, satellite terminal design and on board instrumentation. All these test 
objectives were met. The solar cells worked so well that its transmitters were active for 
six years, (interfering the 108 MHz frequency) well beyond the expected lifetime. The 
scientific results were also a success. As Soviet authorities and even some in the West 
scoffed at the diminutive size of Vanguard 1 and its lack of sophisticated instrumentation, 
it proved to be a very useful tool. Analysis of-the motion of Vanguard 1 established the 
fact that the earth is not spherical but has a bulge, disclosing an unsuspected stress deep 
within the earth. These measurements indicated that there was large-scale convection 
taking place inside the Earth, which supported the emerging theories of continental dnft 
and sea floor spreading. Analysis of the drag exerted by the atmosphere on Vanguard 1 
proved the atmosphere to be far more extensive and variable than previously believed. 
Perturbations in Vanguard 1’s orbit also led to a more refined estimate of the Earths 


oblateness. [Ref. 5] 


2. 1990-1995 


On September 03, 1990 China launched two atmospheric balloons (1990- 
081B/C), weighing 4kg each to measure the magnetosphere. QQW1 was launched into a 
775x804 km orbit at 99 degrees and decayed on March 1, 1991. QQW2 was launched 
into an 833x886 km orbit at 99 degrees and decayed on July 24, 1991. 


Also on two occasions in mid-1990 (February 03, 1994, February 04, 1995), a 
series of spherical objects were released from the shuttle, typically in pairs. These 
Orbital DEbris RAdar Calibration Spheres (ODERACS) were a few centimeters in 
diameter, and were intended to provide calibration for radar echoes. Appendix 2.1 


provides the data on these objects. 
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3. 1995-2000 


A commemorative functioning replica of the onginal Sputnik was deployed from 
the MIR space station (SPUTNIK-40 — aka Sputnik-II, PS-2, RS-17, Sputnik Jr., 1997- 
058C, 24958) during an extravehicular activity (EVA) by Russian cosmonauts Anatoly 
Solovyov and Pavel Vinogradov on November 03, 1997. It was therefore ш а 383х391 
km orbit inclined at 51.6 degrees. It originally arrived via a Progress automated cargo 
rocket on October 9". The satellite was 1/3 scale and weighed only 3kg and was built by 
French students from the l'Aeroclub of France (radio transmitter), staff from the Russian 
Aeronautical Federation (structure), and funded by various sponsors in the space 


industry. It stopped transmitting on the December 29, and decayed on the May 21, 1998. 


The first satellites launched from a submarine were TUBSAT-N (1998-042A) and 
TUBSAT-NI (1998-042B). The 8kg TUBSAT-N and the 3kg TUBSAT-NI (see Fig. 2.2) 
were two nanosatellites launched on July 07, 1998 as a satellite cluster from a Russian 
nuclear powered submarine with a Shtil-1 converted missile in the Barents Sea, and was 
reported to have cost on the order of $100k (US$1998). The satellites were separated in 


orbit via telecommand and were placed into a 400x776 km orbit inclined at 78.9 degrees. 





Figure 2.2 TUBSAT-N (bottom box) / -N1 (top plate) 


The spacecrafts were flat-box shaped with a solar panel on the largest face, and 
each contained three different experimental payloads provided by the Technical 


University of Berlin (TUB): reaction wheel performance, star sensor performance, and 
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store and forward communication. The latter payload consisted of four independent 
communication transceivers for store and forward communication with a baud rate of 
1200 and 2400 baud. Two transceivers worked in the 2m-frequency band, the other two 
in the 70-cm frequency band with FFSK (Fast Frequency Shift Keying) modulation. An 
additional downlink transmitter with 9600 Baud GMSK (Gaussian Minimum Shift 
Keying) modulation was available. The attitude control of TUBSAT-N consisted of two 
magnetic coils, a magnetometer, a reaction wheel and a star sensor. All attitude control 
devices were also developments of the Technical University of Berlin. Electrical power 
was provided by 9 NiCd-battery cells of 5 Ah (SANYO). The battery cells were 
connected serially and provided an unregulated bus voltage from 9 to 13 V. The current 
utilization of the nanosatellites is being used for tracking medium-sized and large 
mammals, stolen cars and to collect data from autonomous buoys for earth environmental 
observation. These buoys are located in the northern Atlantic Ocean near the Canary 


Islands. 


A second sputnik, Sputnik-41 (aka RS-18, 1998-062C), was launched from the 
MIR spacestation into a 313x318 inclined orbit at 51.7 degrees on the November 10, 
1998. A Progress-M40 cargo rocket delivered 1t to MIR on the October 25, 1998. It was 
financed by the Aeroclub de France, and built by French and Russian students. Sputnik- 
41's 200mW transmitter broadcasted pre-recorded voice greetings in French, English, and 
Russian. The spacecraft measured 230mm in diameter, and weighed 4kg. It decayed 


from orbit on the 11th January 1999. 


The third "junior" sputnik, RS-19 (aka Sputnik-Jr. 3rd, 1999-015C), was launched 
from the MIR spacestation by Jean-Pierre Haignere during a spacewalk on the 16th Apml 
1999. A Progress-M41 cargo rocket delivered it to MIR on the 2nd April 1999. Like 
earlier spacecraft in this series, it was supposed to transmit simple messages, however the 
spacecraft was launched "switched off', as advertising messages by a commercial 
company (Swatch) were carried in breach of International Telecommunication Union 


(ITU) regulations regarding amateur bands [Ref 6]. 
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4. 2000 - Present 
a. JAWSAT 


On January 26, 2000 an OSPSLV (Orbital Suborbital Program Space 
Launch Vehicle) carried several microsatellites within a payload adapter called 
JAWSAT, developed jointly by the U.S. Air Force and Weber State University in Utah. 
Those satellites included FalconSat, an experimental satellite built by the U.S. Air Force 
Academy, ASUSat 1, built by students at Arizona State University; the Optical 
Calibration Sphere Experiment, an inflatable 3.5-meter (11.5-foot) balloon built by 
L'Garde for the Air Force Research Laboratory that served as a target for low-powered 
ground-based lasers; and Opal, a Stanford University satellite that, in turn, deployed six 
smaller "picosatellites" built by Santa Clara University, the Aerospace Corporation, and 


ham radio operators. 


In addition to those satellites, JAWSAT included two other payloads that 
remained attached to it after launch. The Plasma Experiment Satellite Test (PEST), 
provided by NASA's Marshall Space Flight Center, studied plasma found at orbital 
altitudes, while Weber State's Attitude Control Platform tested a low-cost three-axis 


stabilization system. The launched satellites are covered below in more detail. 


A tethered set of Picosatellites, DARPA Picosat (OPAL-#1&2), built by 
the Aerospace Corp. at Rockwell and funded by DARPA and UCLA weighing just 0.5kg 
and measuring 100x750x250mm each (see Fig. 2.3), launched from the OPAL 
microsatellite. They were deployed on February 6th, and communications were 
established 24 hours later. Gold strands in the tether were instrumental in radar tracking 
the pair. The spacecraft performed basic tests on MEMS RF switches. The spacecraft 
primary batteries (lithium thionyl chloride) ran out by February 10th. 
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Figure 2.3 Aerospace tethered picosatellites (artist interpretation) 


A picosatellite named after the Greek Goddess of the moon, the 
ARTEMIS Picosat (OPAL-#3&4, Thelma and Louise), weighed just 0.5kg and was built 
in a 5400 hour effort over a period of 10 months by an exclusively female team of 7 
Santa Clara University students. The spacecraft employed a 68HC11 microcontroller, AA 
battery cells and GaAs solar cells. It was designed for one week of operations and was 
not attitude stabilized. It carried a Very Low Frequency (VLF) receiver. The spacecraft 


were deployed on February 11th, and were reported not to be operational. [Ref. 7] 


MASAT, the Miniature Amateur Satellite, and STENSAT, NASA’s 
amateur radio picosatellite built by Goddard Space Flight Center (GSFC), deployed on 
February 12" from the OPAL microsatellite. MASAT (OPAL-#5, JAK) was a 
picosatellite weighting just 0.5kg, and STENSAT (OPAL-#6) weighed in at just 0.25kg. 


Both spacecraft were reported not to be operational. 


ASUSAT-1 was a 5kg nano-satellite designed and built at Arizona State 
University. It carried GPS, a camera, and a radio amateur FM repeater that only operated 
when requested on the uplink. The batteries failed to charge and therefore the satellite 
operated for approximately 15 hours during which telemetry was received, and is no 


longer operational. 
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b. SNAP-I 


SNAP-1 (2000-033C, 26385) is a 6.5kg nanosatellite developed by Surrey 
Satellite Technology Limited (SSTL) and the Surrey Space Centre (SSC) in the UK (see 
Fig. 2.4). The spacecraft was launched as a piggyback ride on June 28, 2000 on a 
COSMOS-3M launcher from Plesetsk, together with the Tsinghua-1 microsatellite 
(China) and the Nadezhda-06 (Russia) primary payload. It was placed into a 684x707 km 
orbit inclined at 98.13 degrees. It carries a remote inspection payload, and an intersatellite 
link to communicate with Tshinghua-1. It carries a Butane propulsion system with a 3m/s 
capability. Further information will be presented on this satellite in Chapter 4 with 
regards to nanosatellite subsystem design. The spacecraft is still reported to be 


operational. 





с. Munin 


Munin (2000-0, 26621A) is а Swedish 6kg satellite (7.5kg including 
separation system) to measure the electron and ion distribution in the aurora ovals and 
was launched into a 698x1800 km orbit. The satellite carries a spectrometer (DINA, 
Detector of Ions and Neutral Atoms), is cubic in shape and measures 200x200x250mm. 
The Swedish Institute of Space Physics (IRF) and the Dept. of Space Physics of Umea 
University (RYP) designed the satellite. Munin was launched from VAFB, together with 
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EO-1 (Earth Observing-1), and SAC-C (‘Satelite de Aplicaciones Cientificas’-C). The 
last contact with Munin was February 12, 2001. 


B. STUDENT SATELLITE PROJECTS 


In recent years, an increased effort to design, build, and operate small satellites 
has taken place in universities and laboratories all over the world. These microsatellites 
and nanosatellites provide numerous flight opportunities for science experiments at a 
fraction of the cost of larger traditional missions. In addition, there has been an 
increasing trend towards international cooperation on space projects. Tomorrow's 
engineers will find themselves working in an international environment for space 
progress to continue. Appendix 2.2 lists the most current educational establishments 


known to be involved in small spacecraft projects. 


1. SSETI - Student Space Exploration & Technology Initiative 


The European Space Agency (ESA) has started this ambitious educational project 
to involve students into the building of satellites: "... The main objective of this initiative 
is to create a network of students, educational institutions and organizations (on the 
Internet) to facilitate the distributed design, construction and launch of (micro)-satellites 
and potentially more complex projects such as a moon-lander [Ref. 8].” The first 


distribution round of the sub-systems resulted in the following distribution: 


e AOCS: Instituto Superior Técnico, Portugal 

e Communication: UNICAL, Cosenza, Italy 

e Ground stations: Instituto Superior Técnico, Portugal 

e Instruments: University of Florence, Italy 

е Lander, Avionics: “| Escola Politecnica Superior at Universitat de Girona, Spain 
e Mechanism: EPFL(Lausanne), Switzerland 


e Mission Analysis: University of Zaragoza, Spain 
е On-Board Data Handling (OBDH): 
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University of Newcastle, UK 


e Power: Euroavia-Napoli, Italy 
e Programmatic: University of Stuttgart, Germany 
e Propulsion: University of Stuttgart, Germany 


e SSETI Infrastructure: Escola Politecnica Superior at Universitat de Girona, Spain 
e Simulation: TU Vienna, Austria 
° Structures/Configuration: 
University of the Basque Country, Spain 
е Structures/Structural calculations: 
Kingston University, UK 
e Thermal: Manchester University, UK 


2 University nanosatellite program 


The Air Force Office of Scientific Research (AFOSR), NASA, and the Defense 
Advanced Research Projects Agency (DARPA) are jointly funding 10 universities with 
grants of $50k/year over two years to design and assemble 10 nanosatellites (~10kg). The 
universities will conduct creative low-cost space experiments to explore the military 
usefulness of nanosatellites in such areas as formation flying, miniature bus technologies, 
enhanced communications, miniaturized sensors, attitude control, distributed satellite 
capabilities and maneuvering. The satellites are planned to launch mid-2002. The 10 
university nanosatellites provide a broad range of technology demonstrations in the areas 
of miniature spacecraft subsystem components and formation flying. There are also 
numerous science measurements and experiments in such areas as GPS scintillation, solar 


wind, magnetic fields, and upper atmosphere ion density. 


The Air Force Research Laboratory (AFRL) is developing a deployment 
structure, securing a launch, and providing such advanced microsatellite hardware as high 
efficiency solar cells and micropropulsion units. NASA Goddard has also teamed with 
the universities to provide approximately $1.2M funding to demonstrate such formation 


flying technologies as advanced crosslink communication and navigation hardware and 
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flight control algonthms. Numerous industry partners are also supporting the universities 


with hardware and design and testing services. 


This program has the potential to provide significant payoff for very modest 
funding by DoD and NASA given the broad university resources being applied and 
support by industry partners. If these flight demonstrations are successful, it is very likely 
government sponsorship can be secured for follow-on launches of nanosatellites built by 


universities and other agencies. 


a. Three Corner SAT 


This project is a joint effort between Arizona State University (ASU), 
University of Colorado at Boulder (CU), and New Mexico State University (NMSU). 
Aptly named Three Corer Sat (3% Sat), the proposed constellation of three identical 
nanosatellites will demonstrate stereo imaging, formation flying/cellular-phone 
communications, and innovative command and data handling. In addition, each 
University in the 3^ Sat constellation has the opportunity to fly an individual unique 


payload should it desire. [Ref. 9] 


b. ION-F 


Utah State University, University of Washington, and Virginia 
Polytechnic Institute are designing and developing a system of three 10-kg spacecraft to 
investigate satellite coordination and management technologies and distributed 
ionospheric measurements. The three will coordinate on satellite design, formation flying 
and management mission development, and science instruments and mission. Advanced 
hardware for distributed space system to be demonstrated includes micro-pulsed plasma 
thrusters (uPPT), gimbaling magnetic attitude control, and an advanced tether system. In 
addition, an Internet based operations center will be designed to enable each university to 
control its satellite from its own remote location. ION-F will focus on mission objectives 


that would benefit TechSat 21 and future Air Force and NASA missions. In addition, 
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industrial (SDL, Primex, Honeywell) support has been identified, including funding for 
students, hardware, and satellite design support. [Ref. 10] 


c. Emerald 


Spacecraft formation flying is an evolving technology with vast scientific, 
military, and commercial potential that ranges from enhanced mission performance to 
radical reductions in operations cost. As part of the TechSat 21 University Nanosatellite 
Program, Stanford University and Santa Clara University are jointly developing 
EMERALD, a low cost, two-satellite mission for validating formation-flying 


technologies. 


Stanford’s SSDL (Space Systems Development Lab) and Santa Clara 
University’s SCREEM (Santa Clara Remote Extreme Environment Mechanisms 
Laboratory) will work as a unified team to develop, construct, test and eventually operate 
the EMERALD spacecraft. The formation flying experiments will be coordinated through 
Stanford’s ARL (Aerospace Robotics Lab) 


The EMERALD Mission 1s divided into three distinct stages that progress 
from a simple single satellite to two free flying satellites in a coarse formation. Using a 
building block experimental strategy, the research payloads first will be characterized in 
isolation. Then, they will be coordinated and combined to permit simple demonstrations 
of fundamental formation flying control functions such as relative position determination 
and position control. At release, the two spacecraft will be stacked together and will 
travel as a single object. This will allow initial checkout, calibration, and some limited 
experimentation. During the second stage of operation, the satellites will separate and a 
simple tether or flexible boom will uncoil, linking the two vehicles. This tethered stage 
will allow full formation flying experimentation including on-orbit relative position 
determination and simple closed loop relative position control using the drag panels and 


microthrusters. 


During the final stage of operation, the tether will be cut in order to permit 


true two-body formation flying for a limited period of time. The tether will have a simple 
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sub-satellite at its midpoint. Upon ground command, the two halves of the sub-satellite 
will separate. Each satellite will retain half of the tether and half of the sub-satellite, 


providing very rough gravity gradient stabilization. [Ref. 11] 


d. UW nanosatellite (Dawgstar) 


The University of Washington Nanosatellite program is a student run 
program to design, build and launch the smallest self-propelled satellite (15kg) to be used 
in a cluster of satellites flown in 2002 as one of the first distributed satellite testbeds in 
space. UW has university partners (Utah State and Virginia Tech), each of which is 
designing a satellite to allow the "cluster" testbed. The focus on the UW program has 
been on distributed control. The sensing will be done using several technologies relating 
to GPS, and a cross-link system designed by APL-JHU. The actuation will be a set of 
eight pulsed plasma thrusters, to be used for both attitude and position control. In 
addition, there are several other important technology developments in this program, 
including: development of uPPTs for both attitude and orbit/formation control, 
development of horizon and sun sensors using small, lightweight, CMOS camera 


technology and distributed ionospheric science. [Ref. 12] 


3. Miscellaneous Projects 


Project Starshine. The Student Tracked Atmospheric Research Satellite is a 
small, optically reflective spherical spacecraft. The Naval Research Laboratory in 
Washington, DC assembled it from eleven hundred sets of aluminum mirror blanks 
machined by Utah technology students and shipped in kits by project officials to schools 
around the world where students polished the blanks. 


The satellite is 48cm (19ın) and was very bright and easily visible to the naked 
eye. Students recorded their observations online while tracking the satellite. The 
satellite’s orbit decayed at a rate proportional to how much the upper atmosphere was 
heated by solar activity, thus monitoring sunspots. Starshine was deployed by NASA 
from a Hitchhiker canister on the Space Shuttle Discovery into a highly inclined low 
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earth orbit on mission STS-96 ın May of 1999, and the satellite re-entered Earth's 
atmosphere on Feb 18, 2000 


Sunsat. Graduate students at the University of Stellenbosch built this South 
African satellite. It was launched as a piggyback payload on a Delta II from Vandenberg 
on February 23, 1999. SUNSAT is a micro satellite, weighing 64 kg, with dimensions of 
45 cm x 45 cm x 60 cm, with an elliptical polar orbit (620 by 850 km). The satellite 
holds amateur radio transponders and several other instruments that allow digital store- 
and-forward capability and a voice ‘parrot’ repeater that is being used primarily for 
educational demonstrations. The unit has two VHF and two UHF transmit-receive 
systems. In addition to amateur radio and school science payloads, Sunsat carries two 
NASA experiments and an experimental push-broom imager capable of taking pictures of 
Earth. The high-resolution imager operates in real time on S-band frequencies. Images 
also can be stored in computer RAM aboard the satellite and then downloaded at lower 
speeds for retrieval by hams and schools. The participants have made a big effort to use 
the project to inspire interest in science and engineering in South African high school 


students. 


CubeSat. The CubeSat is a nano-satellite design from Bob Twiggs of the SSDL 
(Stanford's Space System’s Development Laboratory). The motivation is to develop a 
standard, off-the-shelf-satellite small satellite kit that can be cheaply built, easily adapted 
for different missions, and launched in clusters to lower the per satellite launch cost. The 
CubeSat is about 10cm per side and weighs a kilogram (see next page - Fig. 2.5), 
allowing student groups to be able to build and launch them for around $50k each. 
Eventually, multiple CubeSats will work together in formation to provide the capabilities 
of a single large satellite. The company OSSS (One-Stop-Satellite-Solutions) has been 
formed to commercialize SSDL’s CubeSat. Several college teams are now building 
CubeSats for launch in 2001: CubeSat at Cuesta College Amateur Radio Organization, 
CubeSat at the University of Tokyo, and Cube-sat at Tokyo Institute of Technology. [Ref. 
13] 
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Figure 2.5 Inside a CubeSat (10cm per side) 


PANSAT. The Petite Amateur Navy Satellite (PANSAT) is a small satellite 
designed and built by military officer students, faculty, and staff at the Naval 
Postgraduate School (NPS). The main objective was to support the Space Systems 
Engineering and Space Systems Operations curricula by providing a ‘hands-on’ hardware 
project where exposure to the many facets of a space system development and life cycle 
can be experienced. PANSAT further provides educational training while in orbit through 
a space-based laboratory for officer students at NPS. PANSAT was launched from the 
Shuttle into a low-Earth orbit on the STS-95 Discovery mission as part of the third 
International Extreme Ultraviolet Hitchhiker (IEH-3) experiment. 


The spacecraft itself provides store-and-forward (packet radio) digital 
communications using direct sequence spread spectrum modulation. PANSAT operates 
in the amateur radio 70 cm band with center frequency at 436.5 MHz, a bit rate of 9842 
bits per second and 9 MB of message storage. Amateur radio ground stations will be able 


to utilize PANSAT via a bulletin-board type user interface. 


С. CURRENT SYSTEMS AND INDUSTRY LEADERS 


1: Аегозрасе Согр. 


In an effort sponsored by Defense Advanced Research Projects Agency 


(DARPA), Aerospace Corporation scientists and engineers are collaborating with 
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‚ockwell Science Center and Stanford University to develop miniature low-cost space 
latforms to valıdate microsystems for space applications and advance the development 
f mass-producible, fully functional nanosatellites (see Fig. 2.6). Aerospace supported 
reflight activities for an experiment involving two tiny picosatellites to be deployed 
rom Stanford University's OPAL satellite (see section 2000-present) after its launch by 
пе new Air Force OSPSLV. The picosats were tethered to emulate formation flying 
vithin the range of low-power radios. A picosat mounted on a 50-meter ground antenna 
ormed the third element of a rudimentary constellation. MEMS radio-frequency switch 


rrays were also tested on this mission. 


The mission represents one of several programs for systematic testing and use of 
AEMS in space to be designed and implemented by The Aerospace Corporation. One of 
he MEMS devices, designed and fabricated at The Aerospace Corporation, comprises an 
tray of 19 microthrusters that could be used to onent a nanosatellite. Each of the 19 cells 
epresents a separate thruster like a solid rocket motor on a launch vehicle (see Chapter 


-B for propulsion components). 


The mission closed out February 10, 2000 and was the first of a series designed to 
ralidate MEMS technology and demonstrate the capabilities of mass-produced miniature 
atellites operating in constellations. Another picosat mission was launched July 29, 2000 
iboard the MightySat 2.1 satellite built by the Air Force Research Laboratory (AFRL), 


ind a third more complex "inspector" mission is planned for 2003. 
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Figure 2.6 Artist Conception of orbiting AEROSPACE Nanosatellite 


23 


2. NASA 


The Nanosat Constellation Trailblazer mission is the fourth deep space mission in 
NASA's New Millennium Program. Known as Space Technology 5 (ST-5 for short), the 
mission will attempt to fly three miniature spacecraft high above the Earth. Each of the 
spacecraft is about the size of a ‘big’ birthday cake; 42cm (17 in) across and 20 cm (8 in) 
high, and weighs about 21.5 kilograms (47 pounds). STS5 will attempt to fly three 
separate spacecraft in а constellation, performing coordinated movements, 
communication, and scientific observations as it was a single, larger spacecraft. The 
payoff from this technology demonstration will enable a senes of multi-spacecraft 
missions in the future. Large numbers of small spacecraft are planned in "constellation 
class" missions in the next century to perform in-situ measurements of space weather 
conditions. Space Technology 5 will fly three nanosatellites in a stable major orbit at 


about 200x 35,790 km altitude. 


The spacecraft will be used to test methods for operating a constellation of 
spacecraft as a single system. The mission will also test eight innovative new 
technologies in the magnetosphere. The mission is planned for launch in 2003 as a 
secondary payload on an expendable launch vehicle. The mission is managed by NASA's 


Goddard Space Flight Center in Greenbelt, Maryland and is budgeted at $28 million. 


92 Foreign Universities 


The purpose of the Surrey Nanosatellite Applications Program, by Surrey Satellite 
Technology Ltd (SSTL) located at Surrey University, is to develop practical, modular, 
multi-mission satellite buses within mass ranges of 1-10 kg, and to demonstrate the use of 
miniature electrical and mechanical COTS technologies in space. It also is to provide 
vehicles for the education and training of students in spacecraft engineenng at the 
undergraduate and post-graduate level. The first accomplishment of this was by the 


development of the ultra-low-cost demonstration spacecraft (“SNAP-1’ — see above 2000- 
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present) within a year, by a team comprising undergraduates, Masters and Ph.D. students, 
with the supporting expertise from SSTL and CSER staff. 


Additionally, MicroStructure Technology (MST) at Uppsala University in 
Sweden will launch in a few years the first European nanosatellite. Researchers within 
the Center for Advanced Micro Engineering (AME) are contributing on microsystems 
research hoping to increase the European nanosatellite knowledge base. A nanosatellite 


project has been initiated at AME to promote system oriented research work. 


D. NANOSATELLITES FLYING TOGETHER 


One distinctive concept, which demonstrates the unique capabilities of 
nanotechnology and many of the above small satellite designs, 1s formation flying. The 
next chapter’s thrust will be to present how to use these smaller nanosatellites in 
numerous formations. Many concepts have started to surface on what types of 
formations are better at certain missions, and cluster lifetime orbit analysis simulations 
are quickly narrowing down the required fuel expense and attitude control needed to 


maintain the formation dynamics. 
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Ш. NANOSATELLITE FORMATION DYNAMICS 


A. INTRODUCTION 


In recent years, the innovative idea of distributing the functionality of current 
larger satellites among smaller, cooperative satellites has been sincerely considered for 
various space missions to achieve goals that are not possible or very difficult to 
accomplish with a single satellite. For instance, one possible use for nanosatellites is 
clusters of satellites flying together in formation for high-resolution, synthetic-aperture 
imaging. In this case, groups of nanosatellites are operated cooperatively to act as a 
sparse aperture with an effective dimension larger than can be achieved by a single, 
larger satellite. The system adds flexibility since the formation and therefore aperture size 
and orientation are adjustable on orbit. [Ref. 14] Flying many satellites in formation 
presents flexibility to mission designers given that the individual satellites will be capable 
of repositioning themselves with respect to each other to achieve diverse tasks. By 
accurately computing the preliminary Keplerian orbit elements, the satellites can realize 
the desired close separation and cluster orientation desired to operate the necessary 


missions. 


In particular, NASA and U.S. Air Force have identified multiple spacecraft 
formation flying (MSFF) as an enabling technology for future missions. [Ref 15] Another 
buzzword denoting this division of labor among smaller satellites 1s becoming referred to 
as Distributed Satellite Systems (DSS) in DoD, NASA and the commercial sectors. 
NASA's Mission to Planet Earth and New Millennium programs have acknowledged the 
benefits of satellite formation flying and have incorporated an enhanced formation flying 
experimentation in the Earth Observing System mission. [Ref. 16] The European Space 
Agency has engineered a sophisticated formation for its Cluster mission to study the 
Earth’s magnetosphere, the Orion program is intended for a low cost display of GPS uses 
in formation flying, and the Laser Interferometer Space Antenna (LISA) mission is a 
heliocentric formation-flying mission intended to identify gravity waves. The US Air 


Force Research Laboratory’s TechSat 21 program is a technology demonstration of the 
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‘virtual satellite’ concept. [Ref. 17] The TechSat 21 program was reviewed in Section II 


(above). 


Advances in miniaturization using MEMS technology leads to estimation of 
nanosatellites masses of less than 1 kg. In some instances, simple satellites-on-a-chip 
may have a mass of order 10° kg (picosats). A particular original concept is the use of 
large numbers of such nanosatellites/picosatellites to shape constellations to permit the 
real-time acquisition of distributed information. For example, a spherical constellation of 
nanosatellites has been proposed to provide a real-time, three-dimensional view of the 
magnetosphere. Each nanosatellite can be thought of as a ‘pixel’ of a three-dimensional 
data set. To supply high-quality spatial and temporal resolution for such a mission, large 
numbers of nanosatellites are necessary. Because these ultra-low mass satellites can be 
passive sensors without active orbit control, environmental effects such as air drag will 
form the evolution of the constellation. For MEMS fabricated nanosatellites, with a mass 


of order 0.1 kg or less, constellations could include several thousand elements. 


Other concepts have envisioned large numbers (>10°) of nanosatellites to be 
deployed from a dispenser to provide a continuous planar ring of satellites for 
communication purposes. Such constellations would be fashioned by dispensing 
nanosatellites over a range of orbit radii at the identical inclination to induce differential 
azimuthal motion, therefore forming a homogeneous ring. With a large number of 
satellites in basically random azimuthal locations, communication links grow to be robust 
for on-orbit failures. Furthermore, the constellation 1s short lived because air drag will 
eventually remove all the nanosatellites. Such concepts are attractive for military 
communications, where the constellation is positioned from a single launch vehicle and 
formed to sustain a dedicated remote operation. A related concept requires clusters of 
nanosatellites to be launched into random orbits to piggyback launches on any available 
vehicle. Analysis shows that 400 nanosatellites can provide ~95% coverage of the globe. 
[Ref. 18] The constellation would require the continual deposition of new nanosatellites 


to make up for losses from on-orbit failures and air drag removal. 


Formation flying clusters of satellites also provide for graceful degradation of 


performance during times of satellite failure. If a single large satellite has a system 
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failure, the entire mission is at risk. Ifa single satellite in a cluster fails, the remaining 
satellites in the cluster may maintain mission objectives at a lower performance level. 
The cluster could then be brought back up to mission design specifications or even 


improved with the addition of another inexpensive replacement satellite. 


The following chapter investigates several satellite formation-flying designs, with 
particular emphasis on the projected circular formation. These cluster formations are 
described with regard to their basic orbital equations of motion and then are populated 
through computation from a technique derived from the Aerospace Corporation. [Ref. 19] 
A MATLAB program was written that determined both the mothersat and eight remotes’ 
position and velocity vectors, which were in turn entered into STK. This basic formation 
simulation provides an initial look into the orbital dynamics and a brief look into 


perturbations that affect a LEO cluster with the AV required to maintain that cluster. 


B. INITIAL EQUATIONS OF MOTION OF SATELLITE FORMATIONS 


Satellite formations that are optimal for Interferometric Earth imaging are 


discussed in AIAA paper 98-4379, Optimization of Geosynchronous Satellite Constellations 
for Interferometric Earth Imaging {Ref. 20], and in the following sections four satellite 
formations (in-plane, in-track, circular, projected circular) are obtained for a variety of 








satellite applications. The desired trajectories are carefully designed natural orbits so that 


the energy costs to fly along these trajectories is minimized. 


Satellite formation flying designs can be derived from the linearized equations of 
relative motion for two objects under the influence of a point mass gravitational field. 
These equations are commonly known as Hill’s equations. From Vallado’s text [Ref. 
21], a detailed derivation of Hills’s equation is presented, which take the following form 


for unperturbed motion: 


X-2oy—30°x =0 
y*2ox -0 (3.1) 


¿+02=0 
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Here, x is the radial difference between the two objects, y is the along-track difference, z 
is the cross-track difference (see Fig. 3.1 below), and @ is the mean motion of the 


reference object. 





/ Reference Orbit 
/ 


/ 


Figure 3.1 Hill’s reference frame for satellite relative motion. 


The unperturbed version of Hill’s equations can be solved analytically, with the 


solution being: 


x(t) = a sin(@t) + (3x, + 2%) соѕ(01) + (4х, + 2% 
@ © @ 


y(t) = 2а cos(at) + (6x) + 4); sin(or) - (60x, + 3), )t — 2m y, (3.2) 
4) @ @ 


Ze 
z(t) 2 —-sin(ot) + z, cos(ot) 
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In order to avoid secular, or long-term, growth in the relative motion, Eq. (3.2) needs a 


constraint on its secular term to be set to zero: 


Yo =-2x,0 (3.3) 
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It can then be shown that this constraint results in a displaced orbit with the same energy, 
and thus the same semi-major axis, as the reference orbit through first order in the small 


displacements. By enforcing the constraint, Hill’s equations become: 


al X sin(or) * x, cos(ot) 
0 
| | 2; 
nas 220 cos(@r) — 2x, sin(@t) — Ро У (3.4) 
о о 
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z(t) =—sin(@t) + z, cos(@t) 
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where terms with ‘0° subscripts refer to initial conditions. These equations provide for 


formation flying design. 


In the above equations, it can be seen that the radial (x) and along-track (y) 
components of motion are uncoupled from the cross-track (z) component of the relative 
motion. The motion in the radial/along-track plane can be shown to be an ellipse of fixed 


eccentricity (e = 0.866) with an arbitrary along-track offset. All unperturbed formation 


flying designs must project this elliptical motion into the radial/along-track plane. 


In the lineanzed model, the cross-track (z) component of the relative motion is a 
simple harmonic oscillator. Combining the elliptical motion in the radial/along-track 
direction with the oscillatory motion in the cross-track direction yields the farmly of 


ellipses that describe all formation-flying designs. [Ref. 14] 


Sıx initial conditions must be specified in the solutions to Hill’s equations. These 
initial conditions can be thought of as Cartesian or Keplerian element differences 
between the two satellites and allow for six constraints or design parameters to be placed 
on the formation. One constraint was specified when the secular terms were removed 
from the solution to Hill’s equations (see Eq. (3.3) above). This can be thought of as 
requiring the semimajor axis of the two satellites to be equal. Another constraint 15 Ше 


offset of the elliptical motion in the radial/along-track plane, y,. This leaves four design 


parameters, which describe the size, eccentricity, orientation, and initial location in the 


ellipse of relative motion. These also can be thought of as the size of the ellipse in the 
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radial/along-track plane, the initial location within that ellipse, the amplitude of the 


oscillation in the cross-track plane, and the initial location in the cross-track oscillation. 


1. Basic Remote Motion in the Formation 


The remote formation, as mentioned above, is derived from Hill’s equations. The 
center of the formation is a satellite flying in a circular orbit, defined as ‘mothersat’, 


which then has a moving frame attached to it. Its x and y directions are shown in Fig. 3.2 





Figure 3.2 Moving Frame 


The z direction is pointing upward orthogonal to both the x and y axes. Under 


this coordinate system, the dynamic equations of a remote, “Sat2’, are: 


%-2о)-оҚһад8-" + х) +y +z T} =y, 
j+20x-o’yil-nln +)” +y Haz, (3.5) 


2+ or zin +x) +y F 212 = У, 
Under the moving frame, the desired trajectory of Sat2 in the circular formation satisfies: 


x,(t) = (x, /@)sin(@t) + x, cos(wr) 
y, (t) » 2x, / 9)cos(ot) - 2x, sin(wt) (3.6) 


z,(t) = NS (x, /@)sin(@t) + “Зх, cos(wr) 
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The desired path 1s a nonthrust trajectory for the linearized dynamics of Eq. (3.5). Тһе 


free variables, x, and x, are the initial value of x and the derivative of x. In the inertially 


fixed frame, the path 1s: 


0055 / оа (557) 


All formation designs that follow will be constrained to these equations of 
motion of the remotes (Sat2s). As will be seen, most parts of these equations can be 


canceled out, or simplified. 


2. In-Plane Formation 


The in-plane formation is the simplest of all cluster designs. It consists of a group 
of satellites occupying the same orbital plane and separated by mean anomaly. In Hill's 
equations, setting all initial conditions, except for y,, to zero, represents this formation. 


The solutions are then: 


x(t) 20 
Уб) = у, (3.8) 
20) = О 


where y, represents the amount of in-plane spacing between two satellites. This can be 
related to the mean anomaly separation by: 


Ny > | (3.9) 
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where AM is the mean anomaly separation and a is the semimajor axis of the orbits. 
Again, this formulation being based on Hill’s equations assumes that the orbits are 
circular, are at least nearly so. The primary advantage of the in-plane formation is its 


simplicity in design, deployment and control. 
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3. In-Track Formation 


The in-track cluster design 15 a special case of the in-plane formation. Here the 
satellites all share the same ground track. To do so, the satellites have to occupy slightly 
different orbital planes separated by nght ascension of the ascending node (Q), which 
accounts for the rotation of the earth. The difference in the equations can be seen by the 
addition of a cross-track oscillation that represents the difference in nght ascension of the 


ascending node. The solutions to Hill's equations are then: 





B 
y(t) = У (3.10) 
zin te sin(7) y, cos(@t) 
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where y, represents the amount of in-plane spacing between two satellites, wis the mean 
motion of the orbits, 7 is the inclination of the orbits, and œ, is the rotation rate of the 
Earth. The trailing satellite is behind the lead by some difference in mean anomaly; 
relate this difference in mean anomaly to a difference in time (i.e. At). Then calculate the 
nodal difference AQ such that the second satellite is over the same point as the lead 


satellite at some time in the future (+ + Af). The equations take the following form: 


АМ =2° 
а 
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where AM is the mean anomaly separation and AQ is the nodal separation. The 
attractiveness of the in-track formation is that each satellite in the formation passes over 
the same exact spots on the ground, which is valuable for Earth sensing missions. An 


example screenshot from STK for an in-track formation is shown in Appendix 3.1. 


4. Circular Formation 


The circular formation is one in which satellites maintain a constant distance from 
each other. It can be derived from Hill’s equations analytically or geometrically. The 
analytic approach takes the solutions to Hill’s equations and determines relations between 


the initial conditions given the constraint: 
x+y tz =r? (3.12) 


where r 1s the radius of the formation. The geometric approach takes advantage of the 
fact that the relative motion in the radial/along-track planes (x/y) is fixed in eccentricity. 


From either approach, the following relations are found: 


Di 
Jo == № 
(о 
P 20X; 


(3.13) 
Zo = +\/3х, 
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where the first two conditions set the along-track offset and drift to zero. These 
constraints show that there are two planes in which the circular formation is possible. 
Both intersect the cross-track/along-track plane along the along-track axis but one is 
inclined 30? to that plane and the other is inclined at —30?. If the 30? case is chosen, 


the following solutions to Hill's equations are found: 
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x(t) = =n sin(ot) * x, cos(ot) 
Q 
2% 
y(t) * ——cos(ot) — 2x, sin(ot) (3.14) 
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z(t)= NE) E sin(ot) 4 „Зх, cos(@r) 
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Note that the formation design again has two free parameters: x, and x,. These free 
parameters specify the radius and phasing of one satellite in its circular path about the 
reference satellite. The four other initial conditions were constrained for eliminating 
along-track drift, eliminating the along-track offset, and setting the eccentricity and 


orientation of the ellipse of relative motion. 


If the initial conditions of the solutions to Hill’s equations presented in Eq. (3.14) 


are formulated in terms of the cluster radius and phasing, the following equations arise: 


r 
X 


Зап? (0) +4 


ro „Зчап” (д) 
2 „Зап? (0) + 4 


where г 15 Ше radius of the formation, and @ 1s the phase angle of the satellite measured 


(3.15) 


clockwise from the cross-track(z) direction. 


From Equation (3.14)/(3.15), given a circular-reference orbit and desired cluster 
radius and phasing, it is possible to convert the reference elements from Keplerian to 
Cartesian, after coordinate transformations add the Cartesian differences to get the 
position and velocity vectors of the satellites in the circular clusters, and then convert the 


Cartesian elements back to Keplerian for all the satellites in the cluster. 


The Keplerian element differences for a circular formation are highly dependent 
on the phasing and initial conditions of the reference orbit. For two arbitrary points 1n the 


cluster, there will be differences in inclination (i), right ascension of the ascending node 
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(О), argument of perigee (w), and mean anomaly (M). Generally, satellites in a cluster 
will have the same eccentricity except in the case where the reference orbit is not circular. 
There are conditions, however, where two satellites in the cluster can have the same 


inclination or night ascension of the ascending node. 


The circular formation has two properties which make it attractive: 1) the satellites 
maintain a constant distance from each other, enabling precise measurements for use in 
missions such as interferometry, and 2) unlike the in-plane and in-track formations, the 
circular cluster presents a two dimensional array, increasing dimensional resolution for 
imaging, geolocation potential, and numerous other missions. An example screenshot 


from STK for a circular formation with four remotes is shown in Appendix 3.1. 


5. Projected Circular Formation 


The projected circular formation is very close in design to the circular formation. 
This formation will be covered briefly here in this section, and in more detail with 
another example of equation setup in Section D (satellite formations with constant 


apparent distribution — see below). 


The difference is that the projected circular cluster only maintains a fixed distance 
in the along-track/cross-track (y/z) plane. Another way to describe the projected circular 
formation is to say that when the ellipse of relative motion is projected onto the along- 


track/cross-track plane, it produces a circle. This results in the following constraint: 
y+z’=r’ (3.16) 
where r is the radius of the projected circle. 


Like the circular cluster, the constraint can be applied to the initial conditions to 


the solutions of Hill's equations to produce: 
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Jo eg 
w 
Yo = -20X9 
(3.17) 
20 =+2% 
2, = 42x, 


where the first two conditions set the along-track dnft and offset to zero. There are two 
planes in which the projected circular formation is possible. Both intersect the cross- 
track/along-track plane along the along-track axis but one is inclined 26.565° to that 
plane and the other is inclined at —26.565?. 


If the 26.565? case 1s chosen, the following solutions to Hill's equations are 


found: 


Me) с % sin(@r) + x, cos(@r) 
о 
7 
y(t) = —-cos(@t) — 2x, sin(@r) (3.18) 
2, 


zi zi sin(@t) + 2x, cos(at) 
The formation design again has the same two free parameters, which represent the radius 
and phasing of projected circular path. 


If the initial conditions of the solutions to Hill's equations 1n Eq. (3.17) are 


formulated in terms of the cluster radius and phasing, the following arise: 


ху = > c0s(8) 
19) 


х = = sin(@) 


38 


Like the circular formation design, the Keplerian element differences show up in 
inclination, mght ascension of the ascending node, argument of perigee, and mean 
anomaly, and the variations are dependent on the reference orbit and phasing. The 
primary advantage of the projected circular cluster over the circular cluster is that a fixed 
distance separates the satellites when the formation is projected onto the along- 


track/cross-track plane. This has applications for ground observing missions. 


С. REMOTE CLUSTERS WITH CONSTANT APPARENT DISTRIBUTION 


The work presented in this section continues the work started at the end of the last 
one - the problem of creating a cluster of satellites that would maintain a constant, or 
nearly constant, shape and size when viewed from the Earth. Such a cluster would seem 
to have a constant apparent (as opposed to physical) distribution. In such a formation, all 
satellites would have a clear field of view of the surface, and would remain in sufficiently 
close formation to share their information. Such formations are of interest for large 
distributed-aperture sensing, for example. Another possibility is that of forming clusters 
from many small, inexpensive satellites, each with a particular type of sensor and some 
computing power. Remaining in a close formation would allow the satellites to share 


information and computing ability among themselves. 


For such missions, it might be unnecessary to maintain extremely precise relative 
positions. It might be sufficient to know the relative position accurately, and to remain in 
close enough proximity to allow intersatellite communication. In addition, an effect of 
constant apparent distribution is that the angular dispositions of the satellites relative to 


their ground targets and each other would be constant. 


Since it has been mentioned that one of the desires of nanosatellite designs is to 
reduce the need for station-keeping thrust, orbits in which the natural motion of the 
satellites keeps them in a cluster are required. The need for station keeping would then be 
reduced to eliminating the effects of perturbations on the array (covered below in Section 
E). These effects will be found in the bulk motion of the array and in the relative motions 


of the satellites in the array. For missions in which precise knowledge of position is the 
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primary goal, some cyclic perturbations may be acceptable if they are sufficiently small 


and well understood. 


1. Problem Description 


In what follows, the motion of a satellite will be described with respect both to a 
moving reference point and as seen by an observer on the surface of a sphencal planet. 
The motion of the reference point can be visualized by thinking of it as the position of a 
(possibly hypothetical) reference satellite on a circular orbit of radius R, which will be 
referred to as the reference orbit. This reference point will serve as the origin of a local 
coordinate system. This reference frame was first described by Hill, who derived it in his 


work on the motion of the moon about the Earth. 


Let R be the position of the reference satellite (mothersat), following the reference 
orbit under idealized two-body motion, and r the position of a nearby point (the remote). 
Both R and r are expressed in an inertial frame centered at the planetary center of mass. 
Hill’s coordinate frame moves with the reference point, and rotates such that the x-axis 1s 
aligned with R. The y-axis is tangent to the reference orbit at the reference point, and 
positive in the direction of the orbital motion. The z-axis completes the orthogonal set, as 
was illustrated earlier in Fig. 3.1. Thus, the reference point is also the origin of the 


moving reference frame. 


Hill's frame allows the motion of the actual satellite to be described with respect 
to the reference point. It will be shown that the motion of the satellite in this frame will 
create an orbit; this orbit will be called the relative orbit, as shown in Fig. 3.3 below 
(passing through the plane at some arbitrary angle). The primary interest however, is in 
the relative motion as it appears to a fictional observer on the planetary surface, whose 
position is always on the ray connecting the center of the planet to the reference point. 
This observer represents the point on the surface that is the subject of observation by the 


satellite cluster. 


The concept of an apparent orbit now can be defined. Consider a point on the 


planetary surface that moves with time such that it is always between the center of the 
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planet and the origin of Hill’s frame. The apparent orbit is the motion of the satellite 
relative to the reference point as seen from this point. As this 1s purely a matter of the line 
of sight from the viewer to the satellite, there is no physical meaning to the apparent 
orbit. However, it is helpful to visualize it as the trace left by the intersection of the line 


of sight as it passes through the y - z plane in the Hill’s coordinate frame, as shown in 







Bier 323: 
Relative Orbit 
| 
Reference Point = 
2 
Be Viewer =. 
O 
- 
22 МХ 
Apparent Orbit A 
Reference Orbit — 1. i 
Figure 3.3 Reference, relative, and apparent orbits 


The analysis that follows will be primarily concerned with three angles that 
describe the position of the satellite on the apparent orbit. The elevation angle x, the 
lateral angle A, and the rotary angle $ are shown in Fig. 3.4. The angles are all positive as 


shown, with ó increasing in right-hand rotation about the Hill’s-frame x-axis. 
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Figure 3.4 Geometry description and definitions 


The maximum value of 4 will be referred to as the angular width of the orbit, and 
the maximum value of xas the angular height. Rp will denote the radius of the planet, and 
the value R — Rp will be called the altitude of the reference orbit. This is also the mean 


altitude of the satellite cluster. 


Using this definition, a perfectly circular apparent orbit would result in a constant 
total angle between the line of sight and the line joining the viewpoint to the origin of the 
Hill's frame. The apparent orbit can then be calculated by multiplying the angles by A — 
Rp. The difference in this formal method and the visualization suggested above is 
extremely small for clusters in which the cluster radius is small with respect to its mean 


altitude. 


The symbols a and e will denote the semimajor axis and eccentricity of an orbit. 


The variable @ is the mean motion or average angular motion of a body on an orbit, and 


for an elliptic orbit is given by the relation w= ига? , where yw is the gravitational 
parameter of the central body (the mass of the planet times the universal gravitational 
constant). The true anomaly (v) of a point on an eccentric orbit will be measured from 
periapsis, as shown in Fig. 3.5. In this case, the reference plane will be the plane of the 
circular orbit, and the line of nodes will be taken to be the line of intersection of the two 
orbit planes. The argument of periapsis (w) of the eccentric orbit will be measured from 


this line of nodes. The inclination will be denoted i and will be the angle between the two 
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orbit planes. Finally, the angle 4 2 w-v is the angle from the line of nodes to the radial 


vector. Note that this is defined for circular as well as eccentric orbits. 


Periapsis 






- Reference Orbit 





\ Rising Node 


Dens 35 Orbital elements of the eccentric orbit. 


2: Linearized Approach 


The Clossey — Wiltshire (C-W) equations is an accepted choice for describing the 
motion of a satellite near a circular reference orbit [Ref 21]. They are defined in Hill’s 
reference frame as ın Fig. 3.1, and are generally referred to as Hill’s equations. 


Remember (from Section B) the equations are written: 


Х-2оу-3а/х-0 
ў+20х = 0 (3.1) 


Ž+0 = 
It follows that the out-of-plane motion is given by: 

z(t) = А. соз(@ + $.) (3.20) 
where A; is the magnitude and d, is a constant phase shift. Solving for x and y is only 
slightly more difficult. Following usual practice, integrate y to get: 


y--2ox-k Gay 
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where k is a constant of integration. Substitute this into the x equation to get: 


X--o!x*2ok (802) 


This 1s directly solved to get: 


x(t) = A. cos(@t+¢,)+2k/@ (3.23) 


which is substituted into Eq. (3.1) and integrated to get: 


y(t)=-2A sin(or+¢,)-3kt +e (3.24) 


where c 1s the constant of integration. 


Because a purely cyclic motion is requires with the same period as the reference 
orbit for this method, k = 0; c is simply an offset term and can be given the value zero 


without loss of generality. 


Only the orbit as it appears from the surface of the planet is what is being focused 
upon. In keeping with the assumptions made to linearize the equations of motion, the 
assumption 1s made that the variation in r 1s negligible with respect to the altitude of the 
orbit The appearance of the orbit will then be its projection on the y - z plane. This shows: 

y(t) = 24, sin(ot+0,) 
(3.25) 
Z(t) = A, cos(@t + ¢, ) 


Setting d, = d.results in a perfectly elliptic projection. By properly specifying 
the initial conditions, it is possible to set A, = 2A, and the apparent orbit becomes a circle. 
Setting ø = é + n also results in an elliptic projection, however the direction of the 
projected motion will be reversed. The statement that the apparent orbit is the y-z 
projection assumes that the viewpoint is on the line connecting the center of the planet to 


the center of the relative orbit. 


Because the initial angle ø, in Eqn. (3.25) is arbitrary, it is possible to place as 
many satellites in a circular apparent orbit as desired, each separated from the next by 
some constant angle. This would give a “pinwheel” effect from the planetary surface as 
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the satellites rotate about the center of the cluster. The cluster would spin about its own 


center once during each orbit of the planet. 


The radius of the apparent orbit is also arbitrary. Thus, it is possible to create a 
cluster of an almost arbitrary number of satellites, by creating concentric rings about the 
reference satellite, as in Fig. 3.6. The period of the apparent orbit is always that of the 
reference orbit, so that the entire cluster will maintain its shape relative to the center. 


Thus, separation angles such as $i? in Fig. 3.6 should remain constant. 





Figure 3.6 Apparently circular cluster of satellites. 


As the C—W equations are linearizations, it is to be expected that they do not fully 
describe the relative motion of the satellite. Examining the higher-order terms in the 
orbital equations can make an approximation of their error. This is the subject of the next 


section. 


The height of the apparent orbit can be shown to be independent of the width. 
Thus, this technique can also be used to generate a relative orbit that is very nearly 


circular in actuality, and which appears elliptic when viewed from the planetary surface. 


3. Orbital Mechanics 


The motion of the satellite described by the C—W equations is caused by slight 


differences in the orbital parameters of the satellites. The satellite in the elliptic orbit 
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moves more slowly than the reference when at apoapsis, and more quickly near periapsis. 
It thus seems to lag behind and then to catch up. When combined with inclination, this 


leads to an apparent orbit about the reference satellite. 


Investigating the apparent orbit in terms of the orbital elements is not as straight 
forward as using the linear equations, but this approach allows an investigation of the 
limits of the linearized approach, and an estimate of the corrections due to nonlinearity. 
In this section, the eccentricity and inclination parameters required to produce a circular 


apparent orbit will be investigated. 


In developing the ideas of this section, attention will be confined to orbits in 
which periapsis is 90? from the line of nodes. This ensures the symmetry of the apparent 
orbit about the Hill's-frame z-axis. (Recall that the line of nodes as defined here will be 


the line of intersection of the orbital plane with the reference plane.) 


The description of the apparent orbit will be expressed in terms of the eccentricity 
of the satellite orbit. These terms will be on the order of e; thus, second-order terms in e 
may be considered first-order corrections. In the following analysis, estimates for these 


correction terms will be found. 


a. Apparent Angular Width 


The needed eccentricity is first examined to produce an apparent orbit of a 
particular width. The width will be defined by projecting the line of sight to the satellite 
onto the reference plane. The angle between this projection and the line of sight to the 
reference point is then the angular separation of the satellite and the reference point (this 


Is the angle à in Fig. 3.4). The width of the orbit is then the maximum value of this angle. 


If the inclination of the orbit is zero, the angle is easily computed, taking 
the reference orbit and the eccentric orbit to have equal semimajor axes and thus equal 
periods. The motion of the reference point 1s then equivalent to the mean motion of the 
eccentric orbit. The separation at any point is then the difference between the true 


anomaly v and the mean anomaly M. Assuming the inclination is small enough that its 
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effects are negligible then allows extraction of some information. The expansion as given 


in Vallado [Ref. 21] is: 
у= М = 2esinM +2e” sin2M +... (3.26) 


Finding the maximum, taking the expansion through second order, and taking the 


derivative with respect to M results in: 
4ecos M +5e* cos2M =0 G2) 


For e — 0, this leads to values of M approaching 7/2. Setting M = n/ 2 + ò and dividing 
through by e allows Eq. (3.27) to be rewritten as: 


— 2sinó — 5ecos2ó = 0 (3.28) 


Again, e = 0 gives 6= 0, and for very small e, have 6 is approximately -5e/ 2. The sign of 


ò is less important than the knowledge that ò is of order e. 
Using these results, Eq. (3.26) can now be written: 
v — M 2 2esinó —- $e sin2ó +... (3.29) 


Because 6 1s small, cos ó —1 and sin2ó — 26. Because ò is of order e, the second term on 


the right in the above equation is ~e? and can be ignored. Thus, to second order: 


v-M =2e (3.30) 


Thus, the angular width of the apparent orbit as seen from the center of the central body 
is 4e. From the surface of the planet, the angular spread is wider. Letting A denote the 


angle between the lines of sight to the satellite and to the reference, geometry gives: 


Д = tan”! _rsin(y-M) _ (3.31) 
rcos(v - M)- Е, | 


Expanding the radius of the elliptic orbit in M and e, again writing M as n/ 2 + ò, gives: 
r = RIl-esinö +(e’ /2)(cos2ö -1)+...] | (3.32) 


Using the expansions for v - M and r as expanded, the first order approximation is: 
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(3.33) 


Р 


and that the coefficient of the second-order term is of order 6, making the term of third 


order. This gives a total angular width of the apparent orbit of: 
2А z Ae[ R/(R — R,)] (3.34) 


For a low-altitude orbit, Rp may be nearly as large as R, and thus the angular spread of 
the apparent orbit as seen from the planetary surface might be several times larger than 


the value of v - M. 


b. Correction for Inclination 


To allow for orbital inclination, the expansions used to compute the width 
of the apparent orbit must be extended to include the inclination. Consider the angle 
between the projection of the radial vector onto the reference plane and the projection of 
the eccentricity vector (that is, the vector from the center of the planet to periapsis). 
Letting this angle be denoted v , it shows: 


v= tan” ee) = an ( ВЕ | (3.35) 


rCOSVCOSI COSI 


The width of the apparent orbit is then the difference between this angle 
and the mean anomaly M of the reference orbit. 

This angle will be maximized near v = 2/2. Since the assumption is the 
argument of periapsis is w — 7/2, this corresponds to the descending node. Expanding 
about this point, the small parameter ¢ = v—7/2 is defined, and the angle corresponding 


to V as: 


rsinécosi Е 
tans = —— E Z z= tan НЕ | (3.36) 


rcosé 
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where for brevity use a іп place of cos i. This angle is measured from the line of nodes to 
the projection of the radial vector. The expectation is that £ ~ 0, so that Eq. (3.36) should 


be valid over the area of interest. 


The expansion for £1s: 
e=ae+[(a+a’/3)e? + Ole’) (3.37) 


The angular width is found by introducing the expansion for v into this 
expression. At this point, the assumption is made that the earlier results were “nearly” 
correct, and thus that 6 is in fact of about the same size as e and that i ~ ce for some 


constant c at least through first order. 


Knowing that ò is of the order of e allows the disregard to all terms above 


third order combined in the two variables. Thus: 


=-6=(а-16+(1-а?) аб? /3+2асоз5[1+(1-@?)6?]е 
- [16(1 — a^)ó cos 26 - 5sin2ó6](a/4)e^ —-[32(1 а? )соѕ? 5 (3.38) 
— 3cosó —13cos3ó](a/12) e? +... 


Substituting cos i — cos (ce) for « and expanding the trigonometric terms gives the 


expansion as: 
£—ó6 22e-ó0'e- (54 c^)ó! [2 - (c? + Ye? (3.39) 


through third order combined. Through first order, again the width is 2e, and once more 


there are no second-order terms. 


с. Apparent Vertical Size 


The apparent vertical size of the orbit is determined by the elevation angle 
between the line of sight to the satellite and the reference plane. Relative to the center of 


the planet, this angle can be expressed as: 


к. = ѕіп [sin usin i] (3.40) 
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where u is the angle from the line of nodes to the radial vector. As a result of the 


requirement on the argument of periapsis, the absolute value will be maximized for 


u =. For a circular orbit, it is requires that the maximum value of x, equal the 
maximum value of v — M . Thus, from Eq. (3.30): 
к, -віп (sini)2i-v—M = 2е (3.41) 


to second order in e. 


The vertical angle as seen from the surface of the planet is widened, as is 
the lateral. However, as the maximum vertical spread occurs at the extrema of the radius, 


there is a first-order contnbution of e to r that must be examined. 


When the satellite is at the apses, the apparent elevation angle can be 


written as: 
x = tan” _ RU +e)sini _ (3.42) 
R(l+e)cosi—R, 


where € 1s equal in magnitude to the eccentricity, and is positive at apoapsis and negative 


at periapsis. Expand about ¢= 0 to get: 


ini RR sini 
сі “ое” 6 (3.43) 
Rcosi-R, J (R-R,Y «R'sin i 


Recalling Eq. (3.41) and making small angle approximations, the first term reduces to: 
x =2R-e/(R-R,) (3.44) 


The second term is the first-order correction to x. 


Note that the coefficient of the correction contains R and Rp. An 
assumption is made that the product Re / (R - Rp) is small. This is of concern only when 
Rp approaches R, as for a low-altitude orbit. In this case, it is seen that the correction term 


is of second order in Re / (R - Кр). 


The coefficient of the correction term is negative, showing an expansion 
of the angle at periapsis and a contraction at apoapsis. The effect of this could be to raise 
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the apparent orbit slightly with respect to the reference plane. If in fact the apparent orbit 
is elevated, there should be a corresponding elevation at the points of maximum lateral 


spread. 


To investigate this, a value is required for the position of the satellite on its 
orbit at which the apparent orbit achieves ıt maximum width. Return to Eq. (3.38) in an 
attempt to find this value. If the derivative of this equation is taken, however, no terms 
are found that are of first order combined. Because the value of £- ò, and thus of i, is 
defined only through second order, first-order terms are necessary to obtain any 
information. Numerical investigations show that the elevation of the point of maximum 


width is in fact not linear with e [Ref 21]. 


d. Phase Separation in Apparent Orbit 


When there 1s more than one satellite in apparent orbit about the reference 
point, they will be separated by some angle, constant to the accuracy of the C-W 


equations. In Fig. 3.6, the angle denoted dj; is such a separation angle. 


The separation angle between two satellites in the same apparent orbit is a 
function of the angle between the lines of nodes of their orbits. Nominally, the separation 
angle equals this angle. The separation will vary as a result of the second-order deviations 
of the apparent orbit from perfect circularity; the actual amount of variance will be 


dependent upon the separation. 


Although this makes it impractical to define the change in separation 
angle, the angular rate of the satellite in its apparent orbit can be considered. The satellite 
will sweep its angle faster when it is at periapsis than at apoapsis; the ratio will be the 


same as that for the true anomaly rates: 
6,/6, -[a« o/a- of (3.45) 


where @ is the phase angle in the apparent orbit and the subscripts refer to periapsis and 


apoapsis. 
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An estimate of the total variation in At) from that predicted by the C-W 
equations can be had by finding the value of the true anomaly when the mean anomaly 


equals z/2. This is accomplished through the same expansion used earlier: 
у= М + 2еѕіп М + (5е?/4)ѕіп2М +... у =л[2+ 2е (3.46) 


through second order. The phase angle in the apparent orbit, measured from the vertical, 
can be expressed approximately as: 


rsinusini 


Л 2j 
PM =1/2)=-+ tan | КМ) 


| e 71/2 t sini (3.47) 


Recalling that i « 2e, this implies that the satellite sweeps through 4e less 
apparent angle as M goes from 2А to 5 than during the other half of the orbit. Thus, 


two satellites in orbit about the reference, separated by a 180-degree difference in phase, 


will vary in their relative positions by +4e radians during the orbit. 


е. Geometry of the Orbit 


The apparent orbit described in the previous sections is in a sense the 
projection of the relative orbit onto a plane normal to the radius of the reference orbit (see 
Fig. 3.3/3.4 above). The actual relative orbit is not circular, but elliptical. This follows 


from the results gained in Section C above (Linearized Approach). 


The depth of the motion is clearly the difference between penapsis and 
apoapsis, which by definition is 2Re. As the first-order approximation of the diameter of 
the apparent orbit is 4Re, this implies that the relative motion lies near a plane that is 
angled 60° with respect to the reference plane. As a result of this, the variation in the 
actual distance from the reference point to the satellite will vary from a minimum of the 


apparent orbital radius (when the satellite crosses the reference plane) to a maximum of 


approximately A times this amount (at periapsis and apoapsis) Thus the actual 


maximum distance will be (to first order in e) 5 (R-R,)e. 
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4. Example 


An example to demonstrate numerically the effects of the correction terms on the 
linearized solution is presented. To tie all the chapters together, the nanosatellite design’s 


orbital parameters are used in this example. 


As has been mentioned in the preceding section, the appearance of the relative 
orbit is a matter of the angles between the lines of sight to the satellite and to the 
reference, from a point on the planetary surface along the vector joining the center of the 
planet to the reference. To speak of the radius of the apparent orbit requires a more 
precise definition. For this case, simply multiply the angles describing the relative orbit 


by the altitude of the reference orbit. 


Consider a cluster of Earth satellites, such as in Fig. 3.7, in apparent orbit about a 
central point. The nominal orbital altitude is ~1111 km (600 hr To exaggerate the 
second-order effects, set the radius of the apparent orbit at 25 km. Additionally, assume a 
spherical Earth. The cluster consists of a single ring, with eight satellites equally spaced 


such that their nominal apparent angular separation is 45 deg. 


Nominal 45° 


е о er 


Separation 


Figure 3.7 Ring of eight equally spaced satellites. 


>> 


Using the results derived above, have the angular spread of the desired orbit to be 
A = 25/1111 — 0.0225023 rad. From Eq. (3.34) the equation works out (with R,= 
6378.1363 km, and R=7489.5457 km) that: 


0.0225023 = 2 Re/(R — Rp) => e = 1.67 X 10° rad = .09563 deg 


The inclination angle is the spread angle, so i = 2e = 3.34 x 10° rad = 0.1913 deg. 


The first-order corrections to the vertical spread of the orbit are: 
Ax = 0.110e = 1.84 x 10° = 0.0105 deg 


which when multiplied by the nominal altitude of the orbit produces an apparent shift of 
204.07 m. The shape of the apparent orbit is so close to circular that it would appear so 


to the naked eye. 


D. POPULATING AND MAINTAINING A CLUSTER IN LOW EARTH 
ORBIT (LEO) 


The following analysis was used to generate both cartesian and classical orbital 
elements (Keplerian) that were entered into Satellite Tool Kit? (STK). Both in-plane and 
out-of-plane clusters were developed and varying plane angles of the remotes were 
looked at. Creating a Matlab program (see Appendix 3.2) to calculate the data points, 
some of the final results of the different cluster formations can be found in Appendix 3.3, 


For a cluster constellation such as that in Fig 3.8, the initial orbit elements are 
chosen so that (1) each satellite occupies a node in an arbitrary spatial pattern, and (2) the 
satellites undergo a cyclic motion that allows the formation to persist with relatively little 
AV for maneuvering. Depending on the application, the center position in the cluster may 
be occupied by a mother ship or may be empty. An arbitrary distribution of spacecraft in 
the suborbit plane will maintain its configuration, and the suborbit normal will stay at a 
fixed angle ó relative to the mothersat orbit normal. Large numbers of physically 
disconnected nanosatellites can occupy the suborbit plane and can imitate the 


characteristics of a larger (km-scale) spacecraft. 
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Figure 3.8 Cyclic motion of Subsat Orbit 


I" Populating a Remote Cluster 


Each subsat (remote) undergoes a cyclic motion in the reference frame of a center 
satellite (mothersat). The remote completes one suborbit cycle and returns to its original 
position relative to the mothersat after one revolution about the Earth, as illustrated in 
Fig. 3.8. Therefore all satellites in the cluster must have the same value for the 
semimajor axis (and hence the same orbital period) if perturbing forces are neglected. 
The cluster can be populated based on the idea that each remote is in the 6 o'clock or 12 
o'clock position in the suborbit exactly once per revolution. At this instant the remote 
orbital velocity vector is taken to be parallel to the mothersat velocity vector, but the two 
velocity magnitudes are different. No further assumptions concerning the position and 
velocity of the remote at other points in the orbit are needed to determine the orbital 


elements. 


The initial mothersat orbit is specified by the semimajor axis a, eccentricity ec, 
inclination ij right ascension of the ascending node Q., argument of perigee we, and 


mean anomaly at epoch Moo. Three parameters that specify the cluster geometry (р, и, д) 
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to 


are illustrated in Fig 3.9. The scale factor p determines the overall cluster size of and is 
assigned the same distance value for all remotes (1.e. 25 km), while the dimensionless 
pattern generator 7 can be different for each remote and determines the geometrical 
arrangement (i.e. for circular rings, 7 =1). The angle ô specifies the shape of the suborbit 
and the orientation of the suborbit plane relative to the orbit plane of the mothersat. 
Choosing different values for ó in one cluster will produce a swarm having multiple 


suborbit planes. 
Remote 
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Figure 3.9 Remote Geometry as viewed along Mothersat Velocity vector 


The following paragraphs utilize definitions and equations from Vallado’s 
Astrodynamics textbook. [Ref. 21] Determining the remote orbital elements begins by 
choosing a set of E. values (e.g. 0°, 45°, 90°, etc.) where E. is the eccentric anomaly of 
the mothersat at the location in the orbit when the corresponding remote is at the 6 or 12 
o’clock position (either may be used to populate the cluster). At these ‘population points’ 


the true anomaly, mean anomaly, radius and orbital speed of the mothersat are: 


рено даа eee (E (3.48) 
l-e, 2 


M, =E, —e, sin E, (3.49) 
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етте (3.50) 
COS e 


pu IP 6.50 


where p.=a.(1—e:). From Fig 3.9 the radius and orbital speed of the remote at the 


population point are: 


г, = \/(т їп &)* + (г, + прсо5д) (3.52) 
у, = 42-1) (3.53) 


where the + sign in Eq (3.52) is negative for the 6 o’clock and positive at 12 o’clock. 
Here it is assumed that a, =a,for all remotes, but this choice will be modified later. 


Using the above expressions together with the construction in Fig. 3.9, the radius and 


velocity vectors in the mothersat perifocal coordinates when a remote is at the population 


point are: 
cos@ 
r —r.| sinÓ. 
0 
(3.54) 
-зтд 
У. = > е + с05д 
Pe 0 
(r.+npcosd)cos@. 
r. =| (r.£npcosö)sin®. 
tnpsinó 
(3.55) 


¿Sl 
|! 

E lee 
< 


5 


The magnitudes of these vectors agree with the scalar values in Eqns. (3.50)- 
(3.53). Equation (3.55) expresses the key idea that the remote and mothersat velocity 


vectors are parallel at the population points. 


Knowing that the radius vector and velocity vector of a remote at one point in its 
orbit is sufficient for calculating its orbital elements (as, €s, iss Qs, Ws, Mso). The first step 
is to transform the remote radius and velocity vectors from mothersat perifocal 


coordinates to Earth-centered inertial coordinates as follows: 


7 = К. 7 
(3.56) 
Vi = К. у, 
where R is a 3x3 rotation matrix! defined by: 
R = R(Q,)* R(i,) * Rw,) (3.57) 
and with the three components identified as: 
cos(—-92.)  sin(-2.) 0 
R(2,)=|-sın(-2,) cos(-2,) 0 
0 0 0 
1 0 0 
R(i,)=|0 cos(-i,) sın(-i,) (3.58) 


O =si(=1.) cos(=1.) 


cos(-w) sın(-w) 0 ; 
R(w)-|-sm(-w,) с0$(-м.) 0 
0 0 1 


The second step is to apply the set of equations relating the Earth-centered inertial 


radius and velocity vectors to the classical orbital elements. This procedure yields (ag, és, 


! Matrix found in Matlab program in Appendix 3.2 
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ls, Qs, и;, &) where & is the true anomaly at the population point. The corresponding 


eccentric anomaly and mean anomaly are: 


Е, au | me) (3.59) 
I e. 2 


M,=E.-e,sinE, (3.60) 


The cluster is initialized with the mothersat at perigee, and hence the initial 
mothersat mean anomaly is М№.,=0°. Using Eq. (3.49) for M. at the population point and 


Eq. (3.51)for M; at the same point, the initial remote mean anomaly is M,, — M, — M,, 


which completes the determination of the remote orbital elements. The resulting es, is, 
Q;, w;, Ms. are distributed around the mothersat values, with the difference depending on 


the magnitude of p and on the location on the remote within the cluster. 


The initial subsat orbital elements can be refined to enhance cluster stability in the 
presence of the non-spherical geopotential. This is accomplished by calculating the 
subsat displacements (relative to their original positions) after one orbital period of the 
mothersat. The most significant perturbation (which will be discussed in a later section) 
is due to the second zonal harmonic J2, which appears in the expressions for the mean 


mean-motion w and the secular rates of change of Q and w. 


Propagating the cluster for one revolution in low-Earth orbit under the influence 
of J; reveals that the remote displacements are primarily in the mothersat along-track 
coordinate (y-axis), with cross-track and radial displacements being at least 50 times 
smaller. This suggests that the semimajor axis of each remote can be adjusted slightly to 
cancel the along-track displacement. In the unperturbed solution the semimajor axis was 
the same for all satellites in the cluster, namely a, = a, in Eq. (3.53). The revised 


semimajor axis that compensates for along-track displacements due to J} 1s: 


i mdi NA (3.61) 
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where AY is the along-track displacement after one revolution. Typically the revised 
semimajor axis differs from a, by less than 5m at 1111 km altitude. Of course, adjusting 


a; does not reduce the cross-track and radial displacements caused by J}. 


E. PERTURBATIONS AND AV REQUIREMENTS 


1. STK Perturbation Propagators 


For real analysis the non-spherical geopotential and other perturbing forces 
conspire to disrupt cluster coherency. Data from STK’s numerous perturbation modelers 
propagated the orbits and presented influences based on two-body, J; and Л 
contributions, solar/lunar gravitational effects, solar radiation pressure, and atmospheric 


drag for short term (3-10 days) and long term (1 year) investigation. 


a. Two-Body, J; and J; 


A Two-Body, or Keplerian motion, propagator considers only the force of 
gravity from the Earth, which ıs modeled as a point mass. The two-body propagator uses 
the same basic technique outlined in the two-body equation of motion development. This 
technique assumes the Earth is a perfect sphere and the only force acting on a satellite is 


gravity. This propagator doesn't account for any perturbations. 


J, Perturbation (first-order) and J, Perturbation (second-order) propagators 
account for secular (long-term) variations in the orbit elements due to Earth oblateness. 


These propagators don't model atmospheric drag or solar or lunar gravitational forces. 


J; and Js are zonal harmonic coefficients in an infinite series 
representation of the Earth's gravity field. J? represents the dominant effects of Earth 
oblateness. The even zonal harmonic coefficients of the gravity field are the only 
coefficients that result in secular changes in satellite orbital elements. The J? propagator 
includes the first-order secular effects of the J2 coefficient while the J, propagator 
includes the first- and second-order effects of Jz and the first-order effects of Jz. The J3 


coefficient, which produces long period periodic effects, isn't included in either 
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propagator. J, is approximately 1000 times smaller than J and is a result of Earth 
oblateness. Since the second-order Jz and the first-order J, secular effects are very small, 


there is little difference between the orbits generated by the two propagators. 


b. HPOP 


HPOP is the High Precision Orbit Propagator, and was the main 
propagator for formation analysis for this thesis. HPOP, included as part of STK/PRO, 
can handle circular, elliptical, parabolic and hyperbolic orbits at distances ranging from 
the surface of the Earth to the orbit of the Moon and beyond. As its name implies, it uses 
a powerful propagation technique to incorporate sophisticated orbit perturbation models. 
HPOP uses a Runge-Kutta-Fehlberg integration method of order 7-8 to propagate the 
satellite state in the J2000 reference frame. A variety of high-fidelity models are utilized, 
including the Joint Gravity Model 2 (JGM2; a highly precise model (70 x 70) of the 


Earth's oblateness) with a maximum degree/order of 21. 


The atmospheric density model used the Jacchia-Roberts (similar to 
Jacchia-1971 but uses analytical methods to improve performance) to model drag effects 
on the spacecraft. This model takes into account daily variations in the height of the 
atmosphere due to solar heating among other parameters, and was based on values of 


C, =2.0, daily/average F10.7 of 150, a geomagnetic index of 3.0 and an area/mass ratio 


of 0.020 m?/kg. 


Solar radiation pressure was modeled using С. =1.0 and the same 
area/mass ratio as above. Since sunlight produces a small force on any exposed surface. 
This force varies depending on how reflective the surface is (i.e. a mirrored surface is 
more reflective than a black surface). The solar radiation pressure model in HPOP has 


been updated to be consistent with other commonly used propagators such as GTDS. 


Finally, third body gravity models (solar/lunar gravitational effects) are 
based on U.S. Naval Observatory data and are accurate to within 0.03 arc seconds. 
Additionally, third-body gravitational perturbations are also computed using the DE 405 
planetary ephemeris from JPL for the Sun and Moon. 
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2. Perturbation Effects 


Propagations in the presence of perturbations show the circular formation to be 
highly unstable. The primary factor disrupting the formation design is the earth’s 
oblateness or J» effect. The J contributions to the relative motion are at least an order of 
magnitude larger than the disturbing accelerations including tesseral resonance (for short 
repeat ground track cases), atmospheric drag, solar radiation pressure, and third body 
gravitational effects. Earth oblateness effects are most prevalent in the secular motion of 


the right ascension of the ascending node (£2), argument of perigee (w), and mean 
anomaly (M). 


The cross-track error growth can be attributed to the difference in precession rates 
of the right ascension of the ascending node. Since the two orbital planes have slightly 
different inclinations, the secular J; effect causes the right ascension of the ascending 
node for each orbit to precess at slightly different rates. This results in the orbital planes 


drifting apart and a cross-track error growth. 


An additional contribution to the error growth 1s caused by the rotation of the line 
of apsides of the orbits. This is close to causing the Jz effect to be equal but opposite on 
the mean anomaly and the argument of perigee. If the orbit is not circular (i.e. all the 
remotes’ orbit) the orbit line of the apsides begins to rotate and disrupt the formation 


design. 


Large variations in. eccentricity and argument of perigee that occur in near- 
circular LEO will disrupt the formation, but they can be greatly reduced by placing the 
mothersat into a ‘frozen’ orbit. This choice also simplifies the ephemeris representation, 
because here are no secular or long-period variations in e, and @. The NASA/JPL 
TOPEX mission has successfully demonstrated the long-term stability and maintenance 
of this type of orbit. [Ref. 22] A frozen orbit requires w. = 90° and a particular value of 
e, that is determined by the choice of inclination and altitude. Although only one of the 
satellites in the cluster (typically the mothersat) can achieve frozen conditions, it is 


expected that the total AV of formation keeping will be minimized for this choice. 


62 


Periodic maneuvers must be performed by the mothersat to follow a reference orbit 
having an altitude and frozen eccentricity. Typical the mothersat will have a higher area- 
to-mass ratio then the remotes, and its altitude will therefore decay more quickly. Drag 
compensation is applied so that the mothersat reference orbit matches the average decay 


rate of all the remotes. This minimizes the fuel requirements for each remote. 


For the long-term behavior (~3 years), solar radiation pressure and atmospheric 
drag will disturb the frozen orbit conditions, because the eccentricity and argument of 
perigee are no longer constant. Small oscillations in inclination are common to the 


mothersat and remotes and are not expected to disrupt the formation. 


For the short-term behavior (1-10 days), effects of atmospheric drag, solar 
radiation pressure, and the J to Jg gravitational harmonics disrupt the cluster geometry. 
The position deviations indicate that frequent formation-keeping maneuvers are needed to 
maintain the desired separations among the remotes. Drag is the largest contributor to the 
position deviations at 1111 km altitude, and this is compensated efficiently using 


mothersat maneuvers. The cross-track deviation is caused primarily by differential nodal 


regression: 
-— 5 (3.62) 
3àqo 
лу = 201. (3.63) 


where fm 1s the time between maneuvers, L 1s the in-track tolerance in meters, q is the 
orbit decay rate in m/sec, wis the orbit mean motion. in rad/sec, and AV, is the velocity 
increment of each two-burn transfer. The total AV for the remotes are obtained by 


summing the velocity increments for each maneuver. 


3. Formation Keeping 


The circular and projected circular formations were highly unstable and require 
formation-keeping maneuvers to account for earth oblateness, atmospheric drag, and 


tesseral resonance (for short repeat ground track cycle orbits) effects. 
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Both atmospheric drag and tesseral resonance affect the formation in the along- 
track direction. The along-track error growth induced by these perturbations can be 
controlled via small adjustments in the semimajor axis of the satellites. Based on the 
results of the propagation, these semimajor axis adjustments are sub-centimeter for drag 
effects and on the centimeter level for deep tesseral resonance. From Gauss’ variation of 
parameters (VOP) equations [Ref 23] for Keplerian elements in the normal-tangential 
plane, the change in semimajor axis due to a disturbing acceleration is: 

аа _ 2a°V 


2 3.64 
Р (3.64) 


Ay 


where a is the semimajor axis, V is the magnitude of the velocity, is the gravitational 


constant for the earth, and a; is the acceleration in the direction of the velocity vector. 


If a velocity impulse is assumed, Eq. (3.64) can be rearranged to determine the 


velocity impulse required to produce a desired change in semimajor axis: 


AV, 2 — Aa (3.65) 
Dm 





where the changes in velocity, AV, and semimajor axis, Aa, are assumed small compared 
to the nominal values. From Eq. (3.65), it can be shown that to change the semimajor 
axis by 10cm, a velocity impulse of 0.00478 cm/sec is required for an 1111 km altitude 
near-circular orbit. Thus, the AV, and therefore propellant requirements, to account for 
the differential drag and tesseral resonance effects will be quite small. The frequency of 
the maneuvers will depend on the extent of the drag and resonance effects, formation 
keeping error bounds, and several control related issues, such as the accuracy with which 


these maneuvers are affected. 


The effects of Jz were seen in the cross-track direction with indirect along-track 
contributions. In terms of Keplerian elements, the Jẹ secular effects on right ascension of 
the ascending node, argument of perigee, and mean anomaly are of concern. From Battin 
[Ref. 23], the change in right ascension of the ascending node due to an out of plane 


(cross-track) acceleration is: 
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40 rsin(w+v) 


3.66 
A е 


where r is the satellite radius magnitude, w 1s the argument of perigee, v is the true 
anomaly, h is the angular momentum, i 1s the orbit inclination, and aa 1s the acceleration 
in the direction of the angular momentum vector. If a velocity impulse 1s assumed, Eq. 
(3.66) can be rearranged to determine the velocity impulse required to produce a desired 
change in right ascension of the ascending node: 


Eae: hsin(i) 


y = 49 (3.67) 
r sin(w + v) 


where the changes are assumed small. The differential oblateness effects on nght 
ascension of the ascending node can be derived analytically by taking the partial 
derivative of the governing equation with respect to the inclination (taken from Vallado 


[Ref. 21]): 





2 
Q= кы | А. | wcos(Z) 
Р 


2 

| > (3.68) 
IA E 28 wsın(i)öi 
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Combining Eqs. (3.67)-(3.68), it 1s seen that the amount of velocity 1mpulse 
required to maintain nodal spacing is proportional to the size of the formation and the 


length of the mission: 


ЛО = ОЛ 
2 
m (3.69) 
DEREN = Wr НС 
rsin(w+v)2 p 


where òl is the inclination difference in the formation and Af is the time between 


maneuvers. If the near circular assumption 1s made, Eq. (3.69) simplifies to: 


DIET, 
e BU. (0 ТУСТТ (3.70) 
asin(w+v) 2 
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The total AV requirement is approximately ~7.2m/sec per kilometer of separation per 
year for an 1111-kilometer altitude circular formation. It should be noted that the 
maneuvers must be performed at certain points in the orbit to avoid disturbing the 


inclination. 


The second disruptive influence of the earth’s oblateness on satellite formation 
flying 1s the rotation of the orbit line of apsides. Since the effects of Jz on the argument 


of perigee and mean anomaly are nearly equal and opposite: 


Ber‘ = 
ЛЕ @(4-5sin(i)) 


M, m % vi-e’o@sin’()-2) 
Р 


4 


(3.71) 





The effects of accelerations in the normal direction are also nearly equal and opposite on 


these elements for near circular orbits: 


au = TU + сок | Ain 
dt eV a (3.72) 
dt eaV\a Ж 


In Eq. (3.72), only accelerations in the normal direction (mutually perpendicular to A 
and У are considered since impulses in the tangential (velocity) direction would affect 


semimajor axis maneuvers in the cross track direction do not affect the mean anomaly. 


Since the effects of these maneuvers are nearly equal and opposite on the mean 
anomaly and argument of perigee, the formation keeping analysis can be focused on 
maintaining either one of the elements with maneuvers in the normal direction. This 
assumes that the other element will be maintained by those maneuvers as well. For now, 


the argument of perigee is looked at. 


If impulse maneuvers are assumed, Eq (3.72) can be rewritten in terms of the 


amount of velocity impulse required for a given change in argument of perigee: 


eV 


С Ам (3.73) 


ДИ, = 
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where the eccentricity is considered to be small. The required change in argument of 


perigee can be derived from Eq. (3.71): 


2 
Ду’ = УДГ = IM =) ЛІ 
мие. 


2 
AV =- 4 Ер dc ӘЛІ 
cos(v) 4 p 


(3.74) 





where the amount of AV required is again a function of the length of the mission. For an 
1111 km altitude circular mission, the total Av requirement 15 approximately 2.1m/sec per 
year. It should be noted again that the maneuvers must be performed at certain points in 


the orbit to avoid disturbing the eccentncity. 


4. Station Keeping 


The primary station keeping concerns for all formation flying designs is 
atmospheric drag. For an 1111 km altitude satellite, atmospheric drag could decay the 
orbit's semimajor axis by close to 0.2 kilometers in one year. Of course the effects of 
drag are dependent on the satellite altitude, area to mass, drag coefficient, and the 
atmospheric density, which is a function of the solar cycle. A high-density atmosphere 
near the peak of the solar cycle was used in this analysis so the drag estimates may be 
overly conservative. Station keeping a 0.2 km/year decay in semimajor axis could 


require approximately .097 m/sec of velocity impulse. 


F. OTHER TOPICS 


Another topic specifically important to formation problems is reconfiguration. 
Formations are designed based on the optimization of the various performance metric 
functions for a given mission. These performance metric functions are defined to meet 
the mission requirements. Various primary mission requirements, such as achieving the 
best image quality or gaining the highest probability of detecting moving targets, lead to 


different performance metric functions. As a result, image or interferometry qualities for 
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example, is closely related to the formation baseline and the distribution of the satellites 
in the formation, whereas the performance of moving target indication systems depends 
on the number of satellites and the footprint. To meet multiple mission requirements, it is 
important for the formation controller to have the capability of easy reconfiguration. 
Furthermore, if one satellite has a malfunction in the middle of a mission, the adjustment 
of the satellite distribution to keep the system working, or the replacement of the satellite 
with the malfunction, requires reconfiguration of the formation. Reconfiguration could 
include adjustment of the relative distance between satellites, reassignment of the leader 
of a formation, the changing of the numbers of satellites in the formation, and the 


combination of two formations flying closely. 
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ТУ. SUBSYSTEM DESIGN 


To incorporate a useable and attainable design for potential nanosatellite clusters, 
two satellite subsystems show up as integral to the proposal: 1) the attitude determination 
and control subsystem, and 2) the propulsion subsystem. Although all subsystems of a 
satellite bus are mutually coupled with all other subsystems onboard, these two systems 
allow the close formation flying needed for the numerous cluster missions so far 


mentioned. 


The assumptions offered in the remaining chapter’s component selection are 
derived from given reference payload requirements, the constraint of the previous 
chapter’s calculations of AV required over a course of a year to maintain the cluster 
formation, and to minimize mass, power and physical size. Table 4.1 sums up most 


physical assumptions of the overall spacecraft. 


Characteristic Value 
Satellite shape and composition Cylindrical and homogenous 


Mass 10 kg 


[22 














0.208 m 
0.104 m 


Propulsion & ADCS 
35% (3.5kg) 
Mass Allowance 


Table 4.1 Assumptions made for Satellite Physical Characteristics 
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А. ADCS 


The requirements of ADCS systems are twofold: 
Support formation flight 


e Manage attitude/rates to the extent required to support precision propulsive 


maneuvers 1f available/ required. 


e Maintain pointing within small angles (i.e. + .5° in pitch and roll, + 1° in yaw), 
and maintain rates to within precise control (i.e. +0.1° per second) of the 


desired rates about each axis. 


e Maintain attitude knowledge and spacecraft stability at all times 


Support ground communication 


e Point communication antenna for command uplink and telemetry downlink 


when within sight of ground station 


Attitude Control 1s required for any satellite system and small satellites are no 
exception. Until now most small (i.e. micro-/nanosats) used only very simple and 
rudimentary attitude control systems. This can now change with the development of 
Microelectromechanical Systems (MEMS) and other small consumer electronics devices 
and manufacturing techniques. New control systems are now needed for missions such as 


nanosatellite formation flying missions and space-based interferometry. 


Nanosatellites represents a flexible tool to carry out scientific and technological 
research in space. Nevertheless, obvious limitations in size, mass, onboard available 
power, and costs impose several constraints on the design of nanosatellite subsystems. 
The attitude control subsystem, as one of the more complex subsystems of a satellite, is 
particularly affected by these constraints, especially when the considered applications 
require stringent attitude control. Therefore new solutions in terms of components and 


operating logic need to be investigated to reduce costs, volume, and power requirements. 
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The Attitude Determination and Control Subsystem measures and controls the 
spacecraft’s angular onentation (pointing direction), or, in the case of guidance, 
navigation, and control system, both its onentation and linear velocity (which affects its 
orbit). The simplest spacecraft are either uncontrolled or achieve control by passive 
methods as spinning or interacting with the Earth’s magnetic or gravity fields. These 
may or may not use sensors to measure the attitude or position. More complex systems 
employ controllers to process the spacecraft attitude, and actuators, torquers, or 


propulsion subsystem thrusters to change attitude, velocity or angular momentum. 


1. Magnetic Control 


Because of its low mass and power consumption requirements, magnetic control 
is of particular interest for small satellites, since it is extensively adopted also for active 
attitude damping in gravity gradient stabilized satellites, initial despin, attitude 
acquisition, and procession control. [Ref. 24/25] Several authors have studied and 
documented the attitude control of microsatellites, and a case for nanosatellites can then 
be made from those results. The use of environmental torques, such as those caused by 
gravity, the Earth’s magnetic field, and aerodynamic drag, to control the spacecraft 
attitude by means of passive and/or semi-passive devices allows substantial mass and 
power savings. With regard to this, various solutions have been proposed based on the 
use of gravity gradient booms with eddy current dampers, fluid ring dampers, or soft- 
magnetic damping rods to damp the satellite residual attitude motion. Nevertheless, these 
solutions achieve poor attitude control accuracy (~ 5-10deg). When a finer control is 
required, as in the case of remote sensing applications (~ 0.1 deg), various configurations 
of low-mass, low-power momentum/reaction wheels have been proposed. In this case the 
use of magnetic torquers for momentum dumping, as an alternative to more traditional 
gas jets [Ref. 26], reduces the control system complexity and mass. Components 
available for use in these situations will be presented later on this section (Section A3 


below). 
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а. All-Magnetic Torquer System 


To develop 3-axis attitude control given the very limited power and 
weight availability on a nanosatellite, an all-magnetic torquer system where permanent 
magnets on stepper motors could be used instead of traditional torquer coils. The attitude 
determination would be achieved by a combination of Earth horizon and sun sensors, 
giving three-axis control to approximately two to three degrees. Although this concept 
does not provide fine contro] ability for most remote sensing applications, the idea is to 
progress the knowledge base by getting these nanosatellite systems into space and start 


operational testing and evaluation. 


b. Reaction wheel / Magnetic Torquerod System 


As a last example, the attitude control system of the Italian Scientific 
Microsatellite for Advanced Research and Technology (SMART) was reviewed. Their 
microsatellite attitude control system consisted of three small reaction wheels and three 
magnetic torquers (torquerods). The wheels were used for three-axes attitude control 
during station keeping. The reaction wheel design had been performed using a technique 
that minimized mass and power consumption. As far as the magnetic torquer design was 
concerned, it was primarily driven by the requirements for onboard wheel momentum 
unloading and initial attitude acquisition following the separation from the launcher. 
Nevertheless, the possibility of using the magnetic torquers for attitude contro] during 
station keeping was considered. Their technique was to distribute the torque required for 
attitude control between wheels and magnetic torquers, which minimized the total power 
consumption. The result of this example [Ref. 27] presented the analytical model of the 
proposed technique and demonstrated its effectiveness by numerical simulations. In 


particular, the numerical analysis shows the following: 


e Wheels and torquerods can be simultaneously used to realize a required control 


torque with minimum power consumption. 


e The control torque distribution between torquerods and wheels is governed by the 


torquerod control efficiency, which strongly depends on the Earth’s magnetic 
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field variation along the orbit and on the torquerod design parameter R, / scf 
(where R, is the torquerod winding resistance and scf is the ratio of the torquerod 


magnetic dipole to the torquerod supply current). 


e When more efficient torquerods are used, the control torque portion transferred to 


the torquerods increases so that the total power consumption is minimized. 


e The use of torquerods with low values of the design parameter allows substantial 
power savings with respect to the case of optimal control with reaction wheels 


only. 


e As the required torque percentage assigned to the torquerods increases, the 
attitude control accuracy reduces. Nevertheless, substantial power savings can be 
realized while retaining attitude control accuracy adequate for remote sensing 


applications (0.016 deg for a 40% power savings). 


e The attitude control accuracy could be improved by slightly increasing the 


numerical values of the control gains at the cost of higher power consumption. 


e The proposed contro! technique could be used in addition to wheel momentum 


unloading and attitude acquisition considerations to design the satellite torquerods 


Finally, the proposed technique does not increase the attitude control 
system complexity and mass because the torquerods are also used for initial attitude 


acquisition and wheel momentum unloading. 


2s Propulsion Option 


The attitude control system 1s very closely coupled with several of the other 
systems aboard the spacecraft, and influences or is influenced by every system on the 
spacecraft to some degree. In the remote's attitude control configuration, the most notable 
interaction is normally that with the Propulsion system. One concept that many 
nanosatellite designers are utilizing consists of micro-pulsed-plasma thrusters (uPPTs) to 
provide control of both attitude and translation. The shared nature of the ADCS and 


propulsion actuators also has a profound effect on the formation flight mission, placing 
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limits on the controller design. This design option is explained more in detail under the 


Propulsion section below (Section B1 below). 


3. Components 


Appendix 4.1 shows the overview of ADCS components covered in this section. 


a. Micromechanical gyroscope 


Micromachined solid-state gyroscopes use vibrating mechanical elements 
to sense rotation. They have no rotating parts that require bearings, so they can be easily 
miniaturized. All vibration gyroscopes are based on the transfer of energy between two 
vibration modes of a mechanical structure, caused by Coriolis acceleration [Ref. 28]. The 
highest rotation sensitivity 1s obtained when the drive and sense modes have the same 


resonant frequency. 


Resolution, drift rate, zero-rate output, and scale factor are the most 
important factors that determine the performance of a gyroscope. When a gyroscope 1s 
inertially static, the output signal 1s a random function that 1s the sum of white noise and a 


cyclic noise function of the mechanical resonant frequencies. 


Based on resolution, drift rate, and zero-rate output, the Systron Donner 
QRS-11 micromachined angular rate sensor (see Fig. 4.1) was seen as the leader from 
seven commercially available micromachined gyroscopes to complement magnetometers 
for yaw determination. This gyroscope has a resolution of 0.004°/s, a short-term bias 
stability (100s at constant temperature) of 0.004°/s, and a zero-rate output (DC to 100 
Hz) of 0.01°/s. The gyroscope has a mass of 60g [Ref 29]. 
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Figure 4.1 BEI GYROCHIP™ Model QRS11 Micromachined Angular Rate Sensor 


A close alternative to the QRS-11, would be the Litton G2000 gyro (see 
Fig. 4.2). This component offers a two-axis gyro with a size of only 0.97 inches tall by 
0.76 inches wide and weighs 25 grams [Ref. 30]. This gyro system can be combined 
with a star sensor system to provide attitude position information to provide a pointing 
accuracy of approximately 0.1 degree. The system is still in production and testing and 


has not been space flown. 





Figure 4.2 Litton G2000 gyroscope with electronics 


b. Magnetometers 


Magnetometers are simple, lightweight sensors that measure both the 
direction and magnitude of the Earth’s magnetic field. They are reliable but require 
complex software for interpretation and provide relatively coarse attitude determination 
as compared to horizon, sun, and star sensors. GPS position measurements are used with 
a computer model of the Earth’s magnetic field to approximate the field direction at the 


spacecraft’s current position. Over the course of an orbit, the Earth’s magnetic field 
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direction usually changes rapidly enough with respect to the spacecraft to make 
computation of the field’s time derivative possible. These field variations are large 
enough to enable determination of all three Euler angles using only a three-axis 
magnetometer. The Earth’s magnetic field also varies with time and cannot be calculated 
precisely, so a magnetometer is often used with another sensor such as a sun, horizon or 


star sensor or a gyroscope in order to improve the accuracy. 


The Applied Physics Systems Model 533 miniature three-axis fluxgate 
magnetometer (see Fig. 4.3) can provide direction accuracy to better than 0.1° in a 
laboratory environment. This model is well suited for use in the nanosatellite because of 
its extremely low mass (18g) and its small size [Ref. 31]. Yaw attitude knowledge is 


maintained between magnetometer readings by integration of angular rate measurements. 





Figure 4.3 Model 533: Miniature 3 Axis, Fluxgate Magnetometer 


с: Reaction Wheels 


The reaction wheel is used as an actuator for the attitude control and is one 
of many important key technologies of nanosatellites. Hokkaido Institute of Technology 
(HIT) has developed a small reaction wheel, which is about 150g in weight and 0.015 
Nms in maximum storable angular momentum. The motor, which can be used in a 
vacuum, is 30g and has a torque of 28gcm/2W. The wheel mass balance is adjusted 


precisely and the vibrational level is restricted as low as possible. 
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Figure 4.4 Small reaction wheel developed by HIT 


d. Sensors 
(1) Horizon sensors 


Horizon sensors are infrared devices that detect the contrast between the 
:old of deep space and the heat of the Earth's atmosphere. Horizon sensors can provide 
уќсһ апа roll attitude knowledge for Earth-pointing spacecraft, with an accuracy of 
).1° to 0.25°. For the highest accuracy in low Earth orbit (LEO), it is necessary to correct 


he data for the Earth oblateness and seasonal changes in the apparent horizon. 


Two EDO Barnes Model 13-500 wide-angle miniature solid-state horizon 
iensors (see Fig. 4.5) can be used to provide pitch and roll knowledge to + 1°, including 
ıll errors. These sensors have been space-proven on six missions to date. The sensors 
lave fields of view sufficient to allow pointing at off-nadir angles of up to € 11?, but 
yeak performance is limited to angles less than 9°. Each sensor has a mass of 0.113 kg 
ind is roughly cylindrical with a diameter of about 4.1 cm and a height of about 5.6 cm 
Ref. 32]. 
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Figure 4.5 EDO Barnes Model 13-500 Қаса miniature solid-state horizon 


sensor 


(2) Star Tracker 


Current state-of-the-art commercial star sensors typically weigh 15 
pounds, attain 5 tolO arc-second accuracy, and use roughly 10 watts of power. 
Unfortunately, the current state-of-the-art commercial star sensors do not meet many of 
NASA’s “next-generation” spacecraft and instrument needs. Nor do they satisfy DoD’s 
need for micro/nano-satellite systems. The /ntelligent Star Tracker [Ref. 33], built by 
AFRL, presents a low cost, miniature star tracker for nanosatellite attitude determination 
and navigation. The Intelligent Star Tracker incorporates adaptive optic catadioptric 
telescopes in a single, compact, robust Silicon Carbide housing. The Micro-opto-electro- 
mechanical (MOEMs) micro-mirrors (see Fig. 4.6) are used to compensate for various 
aberrations as well as introduce aberrations such as leveraging off of adaptive optics 
research, the active pixel position sensors enable wide dynamic range and simultaneous 


imaging of faint and bright stars in a single image frame. 
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Figure 4.6 Details of the structure of the micro-mirrors are shown. 


The adaptive optics telescope, using MOEMSs micro-mirrors, facilitates 
extremely accurate tracking, and the ability - when coupled with a star-matching scheme 
based on algebraic coding theory - enables the capability to track at least 5 stars 
simultaneously. Moreover, the massively parallel architecture enables the star tracker to 
operate autonomously without burdening the spacecraft processor and may be used to 
supplement the onboard processor. Because the design utilizes technologies that 
inherently integrate well together and lend themselves to batch processing, estimations 
have that the Intelligent Star Tracker will have a recurring cost less than $100k. In 
addition to low cost, preliminary analysis indicates that the Intelligent Star Tracker will 
have a pointing accuracy exceeding 0.20 arc-sec, NEA better than 0.10 arc-sec, power 


consumption less than 2W and a weight of approximately 200 g. 


е. DGPS 


Formation flying will quickly revolutionize the way science, remote 
sensing and surveillance missions are performed in space, enabling a whole new range of 
applications for small satellites. Currently, there are numerous missions in the planning 
stages involving formation flying of a constellation of micro- or nanosatellites. However, 


to truly achieve the goals of these formation-flying missions, an accurate means of 
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relative ranging, determining time and position measurements,  inter-satellite 


communications, and controlling the formation states is becoming critical. [Ref. 34] 


The need of future formation flying missions to have technology with 
integrated capabilities for communicating, relative ranging, and exchanging precise 
timing among spacecraft within the constellation is fast approaching. AeroAstro Inc. is 
developed a system by integrating a Carner Phase Differential GPS (CDGPS) navigation 
and: attitude sensor with a low power, inexpensive, compact ranging and communications 
system. The result of this integration is a low-cost, robust, secure GPS micro navigation 


and communication system for micro and nanosatellite constellations called Star Ranger. 


The ranging accuracy of Star Ranger is expected to be 1 cm or better, and 
the ultimate goal is set at 3 mm. Using CDGPS, it is expected that the relative attitude 
between spacecraft will be determinable to 0.5? or better. In addition, the overall relative 
position of each spacecraft with respect to each other will also be measurable to less than 


] cm, with a goal of 5 mm. 


B. PROPULSION 


There is an increasing need for on-board propulsion systems for micro- and nano- 
spacecraft. These include upper stage engines to boost spacecraft into final orbits as well 
as very small-scale boosters to launch new 1 kg class payloads for military, commercial 


and scientific applications. 


Decreasing payload sizes will increase demand for smaller, more capable 
platforms, including the ability to maneuver and change orbit; hence the need for small 


propulsion systems. Such propulsive missions could include: 


e Remote inspector to rendezvous and maneuver around a host spacecraft 
e Constellations on the same launch vehicle requiring separation 


e De-orbiting of space junk requiring rendezvous, docking and orbit changing 


In addition to low cost, low mass and short delivery some more specific requirements for 


these propulsion systems include: 
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e Low power consumption 
e Low, controllable thrust 
e High propellant Isp 

e High density Isp 


The simplest spacecraft do not require thrust and hence have no propulsion 
equipment. But most spacecraft need some controlled thrust, so their design includes 
some form of metered propulsion — a propulsion system that can be turned on and off in 


small increments. 


The propulsion system has three objectives. The first objective is to provide the 
remotes with attitude control, which involves rotational disturbance rejection and angular 
positioning. The second objective of the propulsion system is to enable formation flying. 
To maintain a formation with the mother center satellite, and other remotes, the 
propulsion system must have the ability to reject translational disturbances and reposition 
the remotes when the satellite drifts out of the formation. The third objective is to provide 
the AV necessary for the orbital maneuvers throughout the mission as determined by any 
mission plan. These maneuvers could include changing from one formation to another, 
changing orbital parameters, correct velocity errors, maneuver, counter disturbance forces 


(e.g., drag), control attitude during thrusting, and control and correct angular momentum. 


The equipment in the propulsion subsystem includes a propellant supply 
(propellant, tankage, distribution system, pressurant, and propellant controls) and thruster 
or engines. Significant sizing parameters for the subsystem are the total impulse and 


t 


number, orientation, and thrust levels of the thrusters. 


ls Systems 


a. Cold/Hot Gas 


Cold-gas is a more traditional propulsion system. It provides thrust by 
expanding high-pressure gas through a nozzle. Some of the components required for the 


system are tank, tubing, filter, pressure regulator, valves, and thrusters. Currently the 
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GN&C’s Propulsion Branch of the NASA’s Goddard Space Flight Center (GSFC) is 
conducting a broad technology development program for propulsion devices that are 
ideally suited for nanosatellite missions. The goal of the program is to develop 
nanosatellite propulsion systems that can be flight qualified in a few years and flown in 
support of nanosatellite missions. The miniature cold gas thruster technology, the first 
product from the GSFC’s propulsion component technology development program, will 
be flown on the upcoming ST-5 mission in 2003. The ST-5 mission is designed to 
validate various nanosatellite technologies in all major subsystem areas. It is a precursor 
mission to more ambitious nanosatellite missions such as the Magnetospheric 


Constellation mission. [Ref. 35] 


b. MEMS 


Small satellites flying in clusters require periodic “stationkeeping” to keep 
them in place. The required impulse is very small — the goal is not to keep the individual 
satellites in rigid formation, but only to keep them in well-defined orbitals with respect to 
one another. The necessary impulse, therefore, is only the amount needed to overcome 
the difference in drag between the most-affected and the least-affected satellites in the 
cluster. Estimates are that the differential drag can be overcome by providing ~1 mNsec 
(micro-Newton second) to -1 mN sec (milli-Newton second) every 10 to 100 seconds 


throughout each satellite’s mission. 


Currently propulsion technology is developing rapidly towards 
miniaturized systems. Most notable is MEMS technology. The MEMS Mega-pixel 
Micro-thruster Array (see Fig. 4.7) thrusters have very low power and energy thresholds 
for ignition (~10 mWatts, ~100 uJoules), and no moving parts so are expected to be 
highly reliable. [Ref. 36] A single thruster array contains a quarter of a million separate 
thrusters. 
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Figure 4.7 Aerospace MEMS chip compared to Penny 


Marotta Scientific Controls is designing a very tiny 
Microelectromechanical System (MEMS) chip that provides fine attitude adjustments on 
the spacecraft; it uses 8.5 times less power and weighs 2 times less than proven attitude 
control systems. [Ref. 37] Marotta is also in the process of developing a unique 
microthruster (see Fig. 4.8) and electronics driver combination which is capable of low 
power operation (<1 W peak), response time of <5 msec, minimum pulse rate of <1 Hz, 
and has a mass of 50g. The resulting low power component provides an order of 
magnitude reduction in solenoid coil heating when compared to an ordinary continuous 
duty solenoid valve. Aerospace, Primex, Honeywell and AFRL are working separately on 
MEMS based thrusters such as micro-hydrazine. These will be flown on numerous 
commercial and university based nanosatellite systems when the maturity of the 
technology will allow it. The small modular nozzles would allow many options as to 
microthruster size. Although development time will most likely require more than two 


years, the potential for nanosatellites is very high. 
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Figure 4.8 Marotta microthruster compared to Dime 


c. Electrical 


Primex Aerospace Company is working with the University of 
Washington to scale down the power requirements of their full-sized Pulsed Plasma 
Thrusters (see Fig. 4.9). The UW nanosatellite will fly a propulsion system, and will be 


either uPPT”s, or a cold gas system. 





Figure 4.9 Full-sized Pulsed Plasma Thrusters from Primex Aerospace Company 


A typical pulsed plasma thruster consists of two electrodes, a solid 
Teflon® propellant bar, an igniter (spark plug), a feed spring, a power supply, and a 
capacitor (shown in Fig. 4.10). 


The power supply charges the capacitor, which is connected to the two 


electrodes. When a small plasma puff from the spark plug is released between the 
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electrodes, the puff creates a low-resistance arc path, discharging the energy stored in the 
capacitor. This arc ablates a small amount of the Teflon® propellant bar and turns part of 
it into plasma. The current flowing through the arc also creates a magnetic field, and the 


resulting JxB force accelerates the plasma away from the thruster, thus generating thrust. 









reed leniter 5lectrodee 
Spring 4 
(©) Plasma 
P Ж 
oe Electrode 
Capacitor 
Figure 4.10 Basic diagram of a pulsed-plasma thruster 
2. Performance 


Three parameters determine the performance of a propulsion system. These three 
parameters are thrust (T), minimum impulse bit (Jsi), and specific impulse (Zp). Given 
propulsion system hardware, the thrust can be measured experimentally. From these 
experimental data, the average thrust can be calculated. However, for a theoretical 


analysis, the thrust is calculated from Eq. (4.1): 
T=m: Cy (4.1) 


where m is the mass flow rate at the thruster exit and C, is the effective exhaust 


velocity. 


The J;;, 1s the smallest amount of thrust that can be delivered by a thruster. This is 
given in the units of an impulse (force multiplied by time), such as 70 uNs. Although the 
Ij; may be calculated theoretically in some cases, it is usually calculated from the 


experimental data since the minimal thrusting capability depends highly on the 
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propulsion system hardware. The /,;; is calculated using Eqn. (4.2) from the thrust profile 
data of an experiment in which a propulsion system is activated for the shortest duration 


possible. 
Ше (4.2) 


In equation (4.2), t; and ty are the time at the start and the end of the thrusting, 
respectively. For the most spacecraft a low /,;, is desired for fine control of the attitude, 
and especially during dead-band limit cycling. The specific impulse is generally defined 


as shown in equation (4.3): 





I =- | (4.3) 


where g, is the acceleration due to gravity at the surface of the Earth. The specific 


impulse is a measure of the amount of thrust a propulsion system can provide for a given 
propellant mass flow rate. The /,, diminishes for most propulsion systems over their 
lifetime. The specific impulse is essentially a measure of the mass efficiency of a 


propulsion system, so a high Jsp is desired for the propulsion system. 


3. Propulsion System Comparison (Cold gas / uPPT) 


Two main propulsion systems were analyzed and compared. The results of the 


uPPT and cold-gas system performance analyses are summarızed in Table 4.2 below. 
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Propulsion n Ibit Propellant AV Time | Energy 
System 5 ( s Mass per AV | Duration | per AV 
Type ( u (g-s/m) (s2/m) | (J-s/m) 


+ The performance of the uPPT was analyzed assuming a I Hz firing frequency. 





+ The energy per AV requirement for a cold-gas thruster depends on the firing 
mode, pulsed or continuous. 
Table 4.2 Comparison of uPPT and cold-gas propulsion systems (single thruster 


performance). 


Both uPPT and cold-gas thrusters provide enough thrust to compensate for 
maximum translational disturbances, due mostly to the remotes drag. The uPPT system 
has a lower mass due to its lower /;;, and higher 7,,, providing better pointing accuracy 
and more AV for a given propellant mass than the cold-gas system. The fine thrust control 
required for attitude control can be characterized by the long time duration per AV. Thus, 
an impulsive burn, or a short duration per AV burn, which is optimal for most orbital 
maneuvers, 1s of less importance. The peak power and energy consumption are high for 


uPPT systems. 


Most importantly, PPTs do not suffer from propellant leakage. Cold-gas systems 
have a history of failing due to propellant leakage. Compared to the cold-gas system’s 
complexity and problems inherent to a high-pressure system and the limited 
miniaturization capability due to the flow characteristic of gases and liquids, the PPT is 
more attractive. The uPPT has a simple feed system with minimal moving parts, leading 
to higher reliability. Also, PPT system mass can be further reduced by decreasing the 
size of the electronics, the most massive component of Ше иРРТ system. The uPPT 
characteristics are suited to small satellites and thus will become the stepping stone for 
future commercial nanosatellites that will utilize PPTs. 
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EXAMPLE: SNAP-I (SSTL) 


SNAP-1 is a low-cost nanosatellite built by Surrey Satellite Technology Ltd. 
SSTL), it is amongst other objectives a technology demonstrator for 3-axis stabilization 
nd orbit control for a future constellation of small satellites during formation flying. The 
atellite uses a single miniature Y-momentum wheel, 3-axis magnetic torquers and a 
ingle butane gas thruster to ensure a nominal nadir-pointing attitude with full pitch 
ontrol and in-track AV maneuverability. The magnetic torquers do momentum 
jaintenance and nutation damping of the Y-wheel. The primary attitude sensor used, 1s a 
iniature 3-axis fluxgate magnetometer. Precise orbital knowledge was obtained using a 


mall single antenna GPS receiver supported by an on-board orbit estimator. 


A low cost propulsion system was designed and built for the SNAP-1 spacecraft 
ı 7 months from concept to launch site. It utilizes butane stored as a liquid and operating 
1 a cold gas mode. Miniature conventional technology was used for the valves. The 
ropellant was stored in a formed titanium tube, rather than a tank, giving a low cost 
olution. The spacecraft was loaded with 32.6 grams of butane prior to shipping it to 


iunch site. SNAP-1 was successfully launched on 28 June 2000. 


The most obvious feature of the complete propulsion pipework assembly 1s that 
xere was no propellant tank. The propellant was stored in 1.1 meter of coiled titanium 


ibe (see Fig. 4.11), providing 65 cm’ of storage volume. [Ref. 38] 





Figure 4.11 Snap-1 Propellant tube 
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У. CONCLUSION 


А. THE NANOSATELLITE PUSH 


In recent years, an increased effort to design, build, and operate small satellites 
has taken place in universities and laboratories all over the world. These microsatellites 
and nanosatellites provide numerous flight opportunities for science experiments at a 
fraction of the cost of larger traditional missions. This paper has presented an enormous 
thrust toward innovative ways, not only of satellite design, but a shift in thinking about 
satellite’s roles and abilities in the coming years. Nanosatellites are far from being the 
“some distant future’, and indeed are not even the end of the line for this new movement 
іп space commercialization. Picosatellites and femtosatellites are currently being 
designed and built, with space engineering paving the way for many powerful 
opportunities to show how useful and cunning these systems can be to both the 


commercial fields and to the military. 


B. FORMATION DESIGNS 


Imagine satellites numbering in the tens, hundreds or even thousands 
being launched off surface combatant ships, submarines, mobile Army vehicles or even 
Air Force cargo planes. The possibility of throwing a large quantity of ‘sensors’ quickly 
and efficiently into an engaged theater wherever/whenever removes the dependence on 


costly, vulnerable national-asset satellites. 


The preceding chapters looked at the possibility of creating clusters that through 
their natural motion retain a constant shape when viewed from the planet’s surface. It is 
shown that the shape of the apparent circular orbit can be described in terms of 
eccentricity of the orbit, and that terms through second order in e are sufficient to 
describe the motion to a high order of accuracy. Many formation designs have been 
presented to show that utilizing natural orbits allows formation dynamics to help reduce 
fuel requirements for formation-, and station-keeping needs. These natural orbits, 


coupled with robust control laws and precise position knowledge and inter-satellite 
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communication abilities, are key to the growing need and requirement for future NASA, 


DoD and commercial space missions. 


C: PERTURBATION UPKEEP 


This paper has presented the effects of the earth’s oblateness on satellite 
formation flying designs like the circular and projected circular clusters. There are two 
components of motion that must be accounted for: 1) differential changes in the nght 
ascension of the ascending node, and 2) secular changes in the argument of perigee and 
mean anomaly. The cost to maintain relative node spacing is dependent on the size of the 
formation. For a circular cluster of 25 km radius, the cost 1s approximately 7.2m/sec per 
year of velocity impulse. The cost to maintain the argument of perigee, not being a 
function of the cluster size, 1s roughly 2.1 m/sec per year for the circular cluster. If the 
maneuvers cannot be coupled, a cluster like the one presented here could require close to 
10m/sec per year of velocity impulse. Based on the orbit propagations and an assumed 
10% error bounded on the formation, maneuvers would be required every 40 hours or so. 
Other perturbing effects may also require maneuvers but at far less frequency and cost 
than the oblateness induced maneuvers. It should also be pointed out that maneuvers 
may be required in all directions and each satellite in the formation must be able to thrust 
in the along-track, cross-track, and radial directions. The amount of formation keeping 


maneuvers will also vary from satellite to satellite within the formation. 


Station keeping cost for the classes of orbits discussed here is relatively small. 
Only atmospheric drag decay of the semimajor axis is a major concern. Velocity impulse 


to correct for this effect is on the order of .097m/sec per year using high drag conditions. 


Another major issue that impacts the frequency and cost of formation keeping 
maneuvers is that of attitude control. The formation keeping maneuvers discussed above 
require thrusting in the along-track, cross-track and radial directions. The satellite 
dynamics are very sensitive to acceleration in the along-track directions. The largest 
maneuvers are required in the cross-track and radial directions. If the satellite pointing 
has substantial errors when cross-track or radial maneuvers are performed, some 


unwanted acceleration might be applied in the along-track direction with significant 
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consequences. This will undoubtedly drive the frequency and cost of formation-keeping 
higher than what has been calculated here. 

Orbit determination knowledge is another factor that will influence formation 
control. The STK's simulated propagations show that centimeter level differences in 


semimajor axis cause significant along-track error growth over time. The ability to 


determine the orbits to this level of precision is vital for precise and efficient formation 


control. 
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Одегас$ А 


Oderacs В 


Oderacs C 


Oderacs D 


Oderacs E 


Oderacs F 


Oderacs 2A 


Oderacs 2B 


Oderacs 2C 


Oderacs 2D 


Oderacs 2E 


Oderacs 2F 


Launched from 


STS60 9Feb94 


STS60 9Feb94 


STS60 9Feb94 


STS60 9Feb94 


STS60 9Feb94 


STS60 9Feb94 


STS63 4Feb95 


STS63 4Feb95 


STS63 4Feb95 


STS63 4Feb95 


STS63 4Feb95 


STS63 4Feb95 


APPENDIX 2.1 
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International 
Mass (g) Initial orbit @ 51.6 deg 
Designator 
1994-006B 4200 225x463 km @ 51.6 deg 
1994-006C 4200 239x451 km @ 51.6 deg 
1994-006D 500 
1994-006E 500 
1994-0062 5 266x272 km @ 51.6 deg 
1994-006G 5 266х271 Кт @ 51.6 дер 
1995-004С 5000 267x277 km @ 51.6 deg 
1995-004D 4200 323x349 km @ 51.6 deg 
1995-004E 500 
1995-004F ÉS 314x320 km @ 51.6 deg 
1995-004G 1.5 320x334 km @ 51.6 deg 
Not catalogued by 
0.5 
NORAD 
ODERACS Data Table 


Decayed 


10 Feb 94 


10 Apr 94 


3 Mar 95 


6 Feb95 


3 Mar 95 


24 Feb 95 


13 Mar 1995 


29 Sep 1995 


7 Feb 1995 


2 Mar 1995 


27 Feb 1995 


20 Feb 1995 


THIS PAGE INTENTIONALLY LEFT BLANK 
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APPENDIX 2.2 


Educational institutions involved in small satellites 


Increasingly, it 1s becoming possible for educational institutions to be involved in small 
satellites. This is largely due to affordable yet sophisticated technology. This then allows 
satellites to be designed and built within the course duration, or via a series of student 
projects. 


The following list of educational establishements are those known to be involved in small 
spacecraft projects. This can be either in experiments, parts of satellites, or entire 
satellites. Where known, specific spacecraft project names apear in brackets, inclucing a 
year when launched (or expected to be launched). For more details on the projects, please 
refer to the satellite and future mission lists. 


All the links on this page will lead you to the individual institutions 


Europe 


Puniversidad Complutense, Madrid, Spain 
Universidad Politecnica de Madrid, Spain (UPMSAT '95) 
Е „геу Space Centre University of Surrey, Surrey, United Kingdom (UoSAT series '80 
onwards) 
echnical University of Helsinki in Finland (HUTSAT) 
niversity of Leicester, England (CATSAT) 
Brecmnical University of Berlin, Germany (TUBSAT series '90 onwards) 
Buniversity of Bremen, Germany (BREMSAT '94, Abrixas '97) 
Royal Institute of Technology ın Stockholm, Sweden. (KTHSat) 


Binstituto Superior Técnico , Lisbon Portugal. (PoSAT-1) 
nitersitaet Kaiserslautern, Germany (Phase 3D) 


Егер. of Space Physics of Umea University (RYP) (Munin) 











North America 
In the U.S., NASA encourages participation in small missions via the University Explorer 
(UNEX) Programme. 
BUniversity of Illinois at Urbana-Champain. 
BStanford University, , U.S. (SAPHIRE and OPAL) 


BlArizona State University, Tempe, Arizona, U.S. (ASUSAT) 
BUniversity of Arizona, Tuscon, U.S., Students Satellite Program (UASAT) 


Bllowa State Univerity, Ames, lowa, U.S. (ISAT) 
BUniversidad Nacional Autónoma de México, Mexico City, Mexico (UNAMSAT) 
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Мата! Postgraduate School, Monterey, California, U.S. (PANSAT) 

BSan Jose State University, San Jose, California, U.S. (SPARTNIK) 

Bisierra Community College, Rocklin, California, U.S. (MINERVA) 
BUniversity of Alabama, Huntsville, Alabama, U.S. (SEDSAT) 

BUtah State University, Logan, Utah, U.S. (NUSAT '85, Webersat etc CAST) 
Ellohns Hopkins University Applied Physics Laboratory, U.S. 

EiPenn State University, U.S. (SAC-B CUBIC instrument) 

BUniversity of Colorado, U.S. 

Massachusetts Institute of Technology Centre for Space Research, U.S. (HETE) 
EiWeber State Univerity, Ogden, Utah, U.S. (NUSAT '85, WEBERSAT '90 and more) 
BUniversity of New Hampshire, U.S. (CATSAT) 

BBoston University, U.S. (TERRIERS) 

BAnahuac University, Mexico (ANISAT Nano-satellite project) 

BiCarleton University, Canada. (SILA microsatellite) 

BUniversity of Toronto (MOST microsatellite) 
















Africa 
BUniversity of Stellenbosch, South Africa (SUNSAT) 


Asia 
BK orean Advanced Institute of Science /Technology (KAIST), S.Korea (KITSAT series) 
EilMahanakorn University, Thailand (TMSAT) 
National Central University, Taiwan (TUU-Sat) 
BIATSB, Malaysia, (TiungSAT) 
EiTsinghua University, China (Tsinghua-1) 


South America 
EiUniversidad de La Frontera, Temuco, Chile, Involved with the design and construction 
of CESAR-1 
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APPENDIX 3.1 


STK Screen Shots 





In-Track Formation(2 remotes) 
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Circular Formation (4 remotes) 





APPENDIX 3.2 


Matlab Code for STK data input in inertial position/velocity parameters 


¿Cluster Formation 

clear 

ele 

Need Mother parameters 
ac - semimajor axis 
sec - eccentricity 
$ic - inclination 
SRAANC - right ascension of ascending node 
Sargc - argument of perigee 
%Мс - mean anomaly at epoch 
SEC - eccentric anomaly 
$TC - true anomaly 

ҮС - radius 

*$VC - orbital speed 


ac = 7489.137000; %km 


eedeg = 0.00; заед 
ec = ecdeg*pi/180; rad 
icdeg = 63.4; $deg 
Шс = icdeg*(p1/180): trad 
RAANcdeg = 0; =deg 


RAANC = RAANcdeg*pi/180;%rad 
argcdeg = 0; +deg 
argc = argcdeg*pi/180; trad 
Pe = ас (1-ес 2); 


mu = 398600.5; $km 3995-2 


Cluster geometry parameters 


<rho - scale factor; determines overall cluster size/ assigned 
same for all subsats 


<nu - dimensionless pattern generator; determines geometrical 
arrangement 
+delta - angle; specifies shape of suborbit & orientation of 


suborbit plane wrt centersat orbit 
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$rs = subsat radial speed 


%\5 = subsat orbital speed 

заз = subsat semimajor axis 
rho = 100: skm 
deltadeg = 63.4; заеа 
delta = deltadeg*pi/180; “rad 
as = ac; s km 


%First choose eccentric anomaly, Ec, of centersat 
numsatsz8 
for m=0:numsats-1 


Ec= (2*pi/numsats*m) ; sEccentric Anomaly - Equal 
Spacing circular 


%SPattern Generator for Desired Pattern 


naar sCircular Spacing, one loop 


S$Calculte Center-Sat parameters 
Tc=2*atan(sqrt ((1+ec)/(1-ec))*tan(Ec/2)); 
Mc=Ec-ec*sin(Ec) ; 
кс-рей ши ес COS (ME) 7 


ve=sqrt (mu* (2/rc-1/ac) ); 


%+Calculte Sub-Sat parameters (use 12 o'clock insertion position) 
аз=ас; sAssumed for all subsats 
rs=sqrt ((nu*rho*sin (delta) ) ~2+ (rc+nu*rho*cos (delta) )%*2) ; 


vs=sqrt (mu* (2/rs-1/as) ) ; 


Calculate Vectors 
Pcsrc*[cos(TC) во (тео 


Vc=sqrt (mu/pc)*(-sin(Tc) ec+cos(Tc) 0]'; 


Rs=[ (rc+nu*rho*cos (delta) ) *cos (Tc) (rc+nu*rho*cos (delta) )*sin(Tc) 
nu*rho*sin(delta)]'; 


Vs=vs/vc.*Vc; 
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¢Transform subsat radius/velocity vectors from centersat coord to 
Earth-centered inertial coord 


Rxx = 3x3 rotation matrix from RSW to IJK frame 


Rraanc = [cos(-RAANc) sin(-RAANc) 0; -sin(-RAANC) cos(-RAANC) 0; 0 0 
HI 

Rie [100; 0 (cost(i-ic)) (srmn(cwo)); O (-sin(-ic)) (cos(-ic))]; 

Rargc = [cos(-argc) sin(-argc) 0; -sin(-argc) cos(-argc) 0; O O 1]; 


R = Rraanc*Ric*Rargc; 


Rijk=R*Rs 
Vijk=R*Vs 


Rsub(m«1,1:numsats-1)-Rijk'; 
Vsub(m-«1,1:numsats-1)-Vijk'; 


Calculate Orbital Elements 


As should equal Ac for all subsats (Assumed) 


end 


display('SubSat1'); 
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APPENDIX 3.3 


Example orbital Parameters from Matlab to STK coordinate transformation 


Orbital 5=0 а, = 6978.137кт alt=600km Ring radius ~25 km 
Parameters 


oe 534 ЕССЕ БЕ СО ( 

| 0 | 634 | 0003 | 180 | 0.006913 deg | 180 | Тор 
|. 15 | 63.403457deg | 0.003 | 225.012421 49 |0 | 134.987579deg | 

|. 45 | 63.403457deg | 0.003 | 315012421deg | 0 |  44987579deg — 

| 6 | 634 | 0003 | O006913deg | о | 0 Bottom 
| 7.5 | 63.396543deg | 0.003 | 44987579deg | 0 1 315.0124214е9 | 

| 9 | 63.393087deg | 0.003 | 89.975158deg | -0.000001deg |  270.024842deg — 

| 10.5 | 63.396543deg | 0.003 | 134.975 9deg | 0 7 22500242140 — 
NEU | «634 je ООО За NNNM VEDIA we| -0.006913deg |. | 1902 2 


For second outer ring: e=.006 (all other parameters stay the same) 


Orbital 5=634  a,=6978.137km alt=600km Ring radius -25 km 
Parameters 


Eccentricity Mean Anomaly 

— 63.400654 0.00328833 0.211983 deg 359.630879 deg 359.953493 deg Top 
Е 63.06986 0.00325885 270 deg 
г ца 

| 6 | 63.400654 0.00328833 179.788017 deg 0.369121 deg 0.046908 deg Bottom 
ee” lee ee el | ee 
ЩЕ 63.730337 0.00326474 90 deg 
ee A A ==: -— =. 
63.400654 0.00328833 0.211983 deg 359.630879 deg 359.953493 deg 


For second outer ring: e=.006 (all other parameters stay the same) 
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Orbital Е - Ring radius 
ee | ем 6 = 63.4 а. = 7489.137 ктп alt=1111km — pe: 
Position Inclination (°) Eccentricity Arg (°) 


oom 63.400508 re _ ee _ 359.61862 este | шю _ Тор 
p ME 63.058986 0.00300701 m ERE ee Ж 


| ол а жаа 

m 63.400508 0.00300701 359.829235 0.38138 | 180 | Bottom 
Бг = 63.741014 000300707 |  — 270 | 049 | 180 | 

зо аы |. ән ее еи 
63.400508 0.00300701 180.170765 359.61862 | 180 | 


For second outer ring: e=.006 (all other parameters stay the same) 
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APPENDIX 3.4 


RAAN over 1 year(deg) 





Degrees 


| E 9) X ie A 2. 


2001001 2001051 2001101 2001151 2001201 2001251 2001301 2001351 
Date (YYYYMMDD.hhmmss) 


Arg of Perigee over 1 year(deg) - uncorrected 


350 
300 
250 
200 
150 
100 
50 += | 
0 4 | 
2001001 2001051 2001101 2001151 2001201 2001251 2001301 2001351 
Оа{е 
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APPENDIX 4.1 


Sensor Performance Mass Power Unit Units Notes 
р Price 
Micro- 0.004°/s | Use between magnetometer 
Gyroscope | short-term bias and | 0.060 | 0.030 | 5450 М ее 
out regularly using the 
resolution 
magnetometer. 
Magnetometer Es Nominal use is for yaw attitude 
нае 0.018 0.030 $2200 only. Can be used for three-axis 
attitude if necessary. 
Built by 


0113 | 0.030 | $40,000 —— Pitch and roll attitude 


Star Tracker pointing accuracy: | 0.200 < 2.0 $100,000 
0.20arcsec 
NEA: better than 
0.10arcsec 


Comparison of attitude sensors 























Reaction 
Wheel 









Tracker’s processor may be used 







to supplement onboard processor 
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